Quantum phase transitions in a charge-coupled Bose-Fermi Anderson model 
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We study the competition between Kondo physics and dissipation within an Anderson model of 
a magnetic impurity level that hybridizes with a metallic host and is also coupled, via the impurity 
charge, to the displacement of a bosonic bath having a spectral density proportional to u s . As the 
impurity-bath coupling increases from zero, the effective Coulomb interaction between two electrons 
in the impurity level is progressively renormalized from its repulsive bare value until it eventually 
becomes attractive. For weak hybridization, this renormalization in turn produces a crossover from 
a conventional, spin-sector Kondo effect to a charge Kondo effect. At particle-hole symmetry, and 
for sub-Ohmic bath exponents < s < 1, further increase in the impurity-bath coupling results in 
a continuous, zero-temperature transition to a broken-symmetry phase in which the ground-state 
impurity occupancy fid acquires an expectation value {hd)o 7^ 1- The response of the impurity 
occupancy to a locally applied electric potential features the hyperscaling of critical exponents and 
uj/T scaling that are expected at an interacting critical point. The numerical values of the critical 
exponents suggest that the transition lies in the same universality class as that of the sub-Ohmic 
spin-boson model. For the Ohmic case s = 1, the transition is instead of Kosterlitz-Thouless type. 
Away from particle-hole symmetry, the quantum phase transition is replaced by a smooth crossover, 
but signatures of the symmetric quantum critical point remain in the physical properties at elevated 
temperatures and/or frequencies. 
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I. INTRODUCTION 

Quantum impurity models have intrigued physicists for 
more than half a century. 1 In recent years, the focus has 
largely been on models that exhibit quantum phase tran- 
sitions (QPTs). Strictly, these are boundary QPTs at 
which only a subset of system degrees of freedom be- 
comes critical^ Boundary QPTs not only serve as pro- 
totypes for the bulk QPTs encountered (or postulated 
to exist) in many strongly correlated systems*" but in 
certain cases they are amenable to controlled realization 
in quantum-dot setups £ 

Of great current interest are dissipative quantum im- 
purity models that describe a dynamical local degree of 
freedom coupled to one or more bosonic modes repre- 
senting a frictional environment. Experiments on single- 
molecule transistors^ have drawn attention to transport 
through nanodevices featuring electron-phonon interac- 
tions as well as local electron-electron interactions. The 
essential physics of these experiments seems to be cap- 
tured in variants^ ' 9 ' 10 ' 11 ' 12 of the Anderson-Holstein 
model, which augments the Anderson impurity models 
with a Holstein coupling of the impurity occupancy to 
a local (nondispersive) phonon mode. The Anderson- 
Holstein model has been studied since the 1970s in con- 
nection with the phenomenon of mixed valence ; 14 ' 15 ! 16 
and has also been adapted to treat the effect of negative- 
U tunneling centers on superconductivity 17 ' 18 The many 
theoretical approaches that have been applied to these 
models have yielded general agreement that phonons 
serve to reduce the effective Coulomb repulsion between 
electrons in the impurity level, or even to produce an at- 



tractive net electron-electron interaction. Most challeng- 
ing has been the study of simultaneous strong Coulomb 
repulsion and strong electron-phonon coupling. Here, 
the most robust solutions have been provided by an ex- 
tension of the numerical renormalization-group (NRG) 
technique, long established as a reliable tool for tack- 
ling pure-fermionic quantum impurity problemsi 19 ' 20 ' 21 
NRG studies 8 - ' 10 ' 15 have shown that in the one-channel 
Anderson-Holstein model, descriptive of a single molecule 
coupled symmetrically to source and drain leads, increas- 
ing the phonon coupling from zero results in a smooth 
crossover from a conventional Kondo effect, involving 
conduction-band screening of the impurity spin degree 
of freedom, to a predominantly charge Kondo effect in 
which it is the impurity "isospin" or deviation from half- 
filling that is quenched by the conduction band. How- 
ever, even for very strong electron-phonon couplings, the 
ground state remains a many-body Kondo singlet and 
there is no QPT. By contrast, a two-channel model de- 
scribing a single-molecular transistor with a center-of- 
mass vibrational mode exhibits a line of QPTs man- 
ifesting the critical physics of the two-channel Kondo 
models 

An even greater theoretical challenge is posed by quan- 
tum impurities coupled to dispersive bosons. A canonical 
example is the spin-boson model, 22 which describes tun- 
neling within a two-state system coupled to a bosonic 
bath. The model has many proposed applications, in- 
cluding frictional effects on biological and chemical reac- 
tion rates,— cold atoms in a quasi-one-dimensional opti- 
cal trap^ a quantum dot coupled to Luttinger-liquid 
leads^ 5 - and study of entanglement between a qubit 
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and its environment i 26 i In many cases, the dissipa- 
tive bosonic bath can be described by a spectral den- 
sity [formally defined in Eq. ((8]) below] that is propor- 
tional to u) s at low frequencies to. The spin-boson model 
with an Ohmic (s = 1) bath has long been known^ 2 - 
to exhibit a Kosterlitz-Thoulcss QPT between delocal- 
ized and localized phases. The existence of a QPT 
for sub-Ohmic (0 < s < 1) baths was for some years 
the subject of debate. 22 i 28 However, clear evidence for a 
continuous QPT has been provided by the NRG, 2 -^^ 
by perturbative expansion in e = s about the delo- 
calized fixed point and through exact-diagonalization 
calculations £2- 

Theoretical activity has also centered on the Bose- 
Fermi Kondo (BFK) model^ in which an impurity spin- 
i degree of freedom is coupled both to a fermionic band 
of conduction electrons and to one or more bosonic 
baths. BFK models arise in the context of unconven- 
tional heavy-fermion quantum criticality treated within 
extended dynamical mean-field theory (extended DMFT) 
(Ref. l34l ) and have also been proposed to describe quan- 
tum dots coupled either to a noisy environment^ 5 , or to 
ferromagnetic leads^ Studies of BFK models having dif- 
ferent spin rotation symmetry — SU(2), XY, or Ising — 
emplo yin g either expansion^! in e = 1 — s or the NRG 
(Refs.lH and [12) have found continuous QPTs between 
phases exhibiting the Kondo effect and localized phases 
in which impurity spin flips are suppressed by the cou- 
pling to the bosonic bath(s). For exponents < s < 1, 
most evidence suggests that the continuous QPTs of 
the spin-boson and of Ising-anisotropic BFK models are 
equivalent. QPTs outside the spin-boson universality 
class have been identified in dissipative models featuring 
a pseudogap in the electronic density of states4°- 

In this paper, we combine the themes outlined in the 
preceding paragraphs by investigating a charge- coupled 
Bose-Fermi Anderson (BFA) model in which the impu- 
rity not only hybridizes with conduction-band electrons 
but also is coupled, via its electron occupancy, to a bath 
representing acoustic phonons or other bosonic degrees 
of freedom whose dispersion extends to zero energy. The 
model was introduced more than 30 years ago 41 > 42 ' 43 in 
connection with the mixed-valence problem. A spin- 
less version of the model was also discussed in the same 
context. 44 More recently, very similar models have been 
shown to arise as effective impurity problems in the ex- 
tended DMFT for one- and two-band extended Hubbard 
models . 45 ! 46 Hitherto, only limited progress has been 
made toward understanding the physics of such models, 
and we are aware of no study of their possible QPTs. 

Our NRG study of the charge-coupled BFA model with 
bosonic baths characterized by exponents < s < 1 
reveals a crossover with increasing electron-boson (e-b) 
coupling from a spin Kondo effect to a charge Kondo 
effect, very similar to that noted previously in the 
Anderson-Holstein model ^ 10 ' 15 However, under condi- 
tions of strict particle-hole symmetry, further increase in 
the e-b coupling leads to complete suppression of Kondo 



physics at a quantum critical point. Beyond the critical 
e-b coupling lies a localized phase in which charge fluc- 
tuations on the impurity site are frozen. For sub-Ohmic 
baths (0 < s < 1), the QPT is continuous and the nu- 
merical values of the critical exponents describing the 
response of the impurity charge to a locally applied elec- 
tric potential demonstrate that the transition belongs to 
the same universality class as that of the spin-boson and 
Ising BFK models. For Ohmic baths (corresponding to 
s = 1), the QPT is found to be of Kosterlitz-Thoulcss 
type. Particle-hole asymmetry acts in a manner analo- 
gous to a magnetic field at a conventional ferromagnetic 
ordering transition, smearing the discontinuous change 
in the ground-state as a function of e-b coupling into a 
smooth crossover. Signatures of the symmetric quantum 
critical point remain in the physical properties at elevated 
temperatures and/or frequencies. 

It is important to note that questions have been raised 
as to whether or not the NRG method reliably captures 
the quantum critical behavior of the spin-boson and Ising 
BFK models for bath exponents < s < i. It is 
a standard belief 3 -^ that the low-energy behavior near 
a quantum phase transition in d spatial dimensions is 
equivalent to that of a classical transition in d + z di- 
mensions, where z is the dynamical exponent. In the 
case of the spin-boson and Ising BFK models, for which 
d = and z = 1, the corresponding classical system 
is a one-dimensional Ising chain with long-ranged inter- 
actions that decay for large separations r like r~( 1+s \ 
The Ising chain is known to possess an interacting crit- 
ical point for | < s < 1, but to exhibit a mean-field 
transition 4 ^ for < s < i. By contrast, NRG studies of 
the spin-boson 3 -^ and Ising BFK (Refs. [H and[H) mod- 
els have found non-mean-field behavior extending over 
the entire range < s < 1, leading to a claim of break- 
down of the quantum-to-classical mapping. 33 * This claim 
has recently been contradicted by continuous-time Monte 
Carlo^ - and exact diagonalization 3 ^ studies. Debate is 
ongoing concerning the interpretation of these various 
results . 50 ! 51 The eventual resolution of this debate may 
determine the validity of the small subset of our NRG re- 
sults that concerns the critical exponents of the charge- 
coupled BFA model with bath exponents < s < A. 
There is every reason to believe that the remaining re- 
sults are physically sound. 

The rest of this paper is organized as follows. Section 
|H1 introduces the charge-coupled BFA Hamiltonian and 
summarizes the NRG method used to solve the model. 
Section IIIII contains a preliminary analysis of the model, 
focusing on the bosonic renormalization of the effective 
electron-electron interaction within the impurity level. 
Numerical results for the symmetric model with sub- 
Ohmic (0 < s < 1) dissipation are presented and inter- 
preted in Sec. IIVI Section [V] treats the symmetric model 
with Ohmic (s = 1) dissipation. Section [VTl discusses the 
effects of particle-hole asymmetry. The paper's conclu- 
sions are presented in Sec. IVIII 
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II. MODEL AND SOLUTION METHOD 



A. Charge-coupled Bose-Fermi Anderson 
Hamiltonian and related models 



In this work, we investigate the charge-coupled Bose- 
Fermi Anderson model described by the Hamiltonian 

-HCCBFA = #imp+^band + ^bath + ^imp-band + -ffimp-bath, 

(1) 

where 



-ffimp = tdfld + Ufld^fldl, 

H band =^2^C ka C ka , 
k,cr 

-ffbath = ^2 w q a q a q , 



-^imp-band 
-^imp-bath 



fNk 



j= {hd - 1) A qK + a -q) 



(2) 
(3) 

(4) 
(5) 
(6) 



Here. 



d a annihilates an electron of spin z component 
a = ±^ (or a = t, 4) and energy < in the impu- 
rity level, h da = d^dl, h d = n dT + h dl , and U > is 
the Coulomb repulsion between two electrons in the im- 
purity level. Vk is the hybridization between the impu- 
rity and a conduction-band state of energy £k annihilated 
by fermionic operator Ck CT , and A q characterizes the cou- 
pling of the impurity occupancy to bosons in an oscillator 
state of energy w q annihilated by operator a q . N/. is the 
number of unit cells in the host metal and, hence, the 
number of inequivalent k values. Correspondingly, N q is 
the number of oscillators in the bath, and the number 
of distinct values of q. Without loss of generality, we 
take Vk and A q to be real and non-negative. Throughout 
the paper, we drop all factors of the reduced Planck con- 
stant H, Boltzmann's constant ks , the impurity magnetic 
moment gfJ-B, and the electronic charge e. 

To focus on the most interesting physics of the model, 
we assume a constant hybridization Vk = V and a flat 
conduction-band density of states (per unit cell, per spin- 
z orientation) 



\po = (2D)- 1 for|e|<£ 
| I otherwise, 

(7) 

defining the hybridization width T = irp$V 2 . The 
bosonic bath is completely specified by its spectral den- 
sity, which we take to have the pure power-law form 



{Kq\) 2 9}- s uj s for < u < ft 







otherwise, 



(8) 



characterized by an upper cutoff Q, an exponent s that 
must satisfy s > —1 to ensure normalizability, and a 
dimensionless prefactor A"oA. In this paper, we present 
results only for the case f2 = D in which the bath and 
band share a common cutoff. We also adopt the conven- 
tion that Kq is held constant while one varies A, which 
we term the electron-boson (e-b) coupling. It should be 
emphasized, though, that the key features of the model 
are a nonvanishing Fermi-level density of states p(0) > 
and the asymptotic behavior B(ui) oc lo s for lu — > 0. Re- 
laxing any or all of the remaining assumptions laid out in 
this paragraph will not alter the essential physics of the 
model, although it may affect nonuniversal properties, 
such as the locations of phase boundaries. 

For many purposes, it is convenient to rewrite^ the 
impurity part of the Hamiltonian (dropping a constant 
term e d ) 



H imp = S d (h d - 1) + — (ha - 1 



(9) 



where S d = e d + U/2. Most of the results presented below 
were obtained for the symmetric model characterized by 
e d = —U/2 or 8 d — 0, for which the impurity states 
n d = and n d = 2 are degenerate in energy. Section IVT1 
addresses the behavior of the asymmetric model. 

In any realization of -ffccBFA involving coupling of 
acoustic phonons to a magnetic impurity or a quantum 
dot, the value of the bath exponent s will depend on 
the precise interaction mechanism. However, phase space 
considerations suggest that any such system will lie in 
the super-Ohmic regime s > 1. Models closely related to 
-Hccbfa have also been considered in the context of ex- 
tended DMFT,-^ a technique for systematically incor- 
porating some of the spatial correlations that are omitted 
from the conventional DMFT of lattice fermions. 48 Ex- 
tended DMFT maps the lattice problem onto a quantum 
impurity problem in which a central site interacts with 
both a fermionic band and one or more bosonic baths, 
the latter representing fluctuating effective fields due to 
interactions between different lattice sites. The charge- 
coupled BFA model serves as the mapped impurity prob- 
lem for various extended Hubbard models with nonlocal 
density-density interactionsi 45 ' 46 In these settings, the ef- 
fective bath exponent s is not known a priori, but is de- 
termined through self-consistency conditions that ensure 
that the central site is representative of the lattice as a 
whole. The extended DMFT treatment of other lattice 
models 34 gives rise to exponents < s < 1, and we ex- 
pect this also to be the case for the extended Hubbard 
models. 

At the Hartree-Fock levelj 4 ^ the impurity properties of 
Hamiltonian {1} are identical to those of the Anderson- 
Holstein Hamiltonian, 



i?AH = H a + LUaa^a + \ {h d - l)(a + a)), 



(10) 



which augments the well-studied Anderson impurity 
models 3 - 



H A = Hn 



Hi 



imp-band; 



(ii) 
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with a Holstein coupling of the impurity charge to a sin- 
gle phonon mode of energy uj . At several points in the 
sections that follow, we compare and contrast our results 
for -ffccBFA with those obtained previously for Hah- 



B. Numerical renormalization-group method 

We solve the charge-coupled BFA model using the 
NRG method ) 19 ' 20 ! 21 as recently extended to treat mod- 
els involving both dispersive bosons and dispersive 
fermions . 38 ' 39 The full range of conduction-band ener- 
gies —D<e<D (bosonic-bath energies < u> < fl) 
is divided into a set of logarithmic intervals bounded by 
e = ±DA~ k (oj = QA~ fe ) for k = 0, 1, 2, where A > 1 
is the Wilson discretization parameter. The continuum 
of states within each interval is replaced by a single 
state, namely, the particular linear combination of band 
(bath) states within the interval that enters -ffimp-band 
(-ffimp-bath)- The discretized model is then transformed 
into a tight-binding form involving two sets or orthonor- 
malized operators: (i) f n(J (n = 0, 1, 2, . . .) constructed 
as linear combinations of all Cko- having \ey\ < DA~ n ; 
and (ii) b m (m = 0, 1, 2, . . .) mixing all a q such that 
< w q < f2A~ m . This procedure maps the last four 
parts of Hamiltonian |T]) to 

oo 

-^band = t \ e nfna-fnt7 "b T n {fncrfn- 1 ,a "b fn - l.afna) 1 i 

(12) 

oo 
m=0 

(13) 

AT P G ba„d = \j^-Uld a + 4/o.) , (14) 

F^ p G bath = nKo \ (n d 1)(6 + bt). (15) 
\/n(s + 1) 

Here, t = to — 0, while the remaining coefficients e n , 
r„, e m , and t m , which include all information about the 
conduction-band density of states p(e) and the bosonic 
spectral density B(u>), are calculated via Lanczos recur- 
sion relations. 39 For a particle-hole-symmetric density of 
states such as that in Eq. |J7J, e„ = for all n. 

The coefficients r„ in Eq. (fT2"]) vary for large n as 
DA~ n / 2 , while e m and t m entering Eq. (|13|> vary for large 
m as J7A~ m . Therefore, the problem can be solved iter- 
atively by diagonalization of a sequence of Hamiltonians 
Hm (N = 0, 1, 2, . . .) describing tight-binding chains 
of increasing length. At iteration N > 0, Eq. (fT2")) is 
restricted to < n < N, while Eq. p^]) is limited to 
< m < M(N). The spirit of the NRG is to treat 
fermions and bosons of the same energy scale at the same 
iteration. Since the bosonic coefficients decay with site 
index twice as fast as the fermionic coefficients, after a 



few iterations the iterative procedure requires extension 
of the bosonic chain only for every second site added to 
the fermionic chain. In this work, we have chosen for 
simplicity to work with a single high-energy cutoff scale 
D = Q,. It is then convenient to add to the bosonic chain 
at every even-numbered iteration, so that the highest- 
numbered bosonic site is M(N) = [N/2\, where [^J is 
the greatest integer less than or equal to x. 

The NRG method relies on two additional approxima- 
tions. Even for pure-fermionic problems, it is not feasible 
to keep track of all the eigenstates because the dimension 
of the Fock space increases rapidly as we add sites to the 
chains. Therefore, only the lowest lying N s many-particle 
states can be retained after each iteration. The presence 
of bosons adds the further complication that the Fock 
space is infinite-dimensional even for a single-site chain, 
making it necessary to restrict the maximum number of 
bosons per chain site to a finite number Nf,. Provided 
that N s and Nb are chosen to be sufficiently large (as 
discussed in Sec. lIVAl) . the NRG solution at iteration N 
provides a good account of the impurity contribution to 
physical properties at temperatures T and frequencies u 
of order DA~ N / 2 . 

Hamiltonian {T]) commutes with the total spin-z oper- 
ator 

Sz = \{n d] - hdi) + ^(A^t ~ /»j/rj' ( 16 ) 
n 

the total spin-raising operator 

S + = d\d l+ J2fi,f nl ^(S-)\ (17) 

n 

and the total "charge" operator 

Q = n d 1 + £(4 T /»t + flf nl - 1), (18) 

n 

which measures the deviation from half-filling of the total 
electron number. One can interpret 

h=\Q, i+ = -d\d\ + $>ir/y^ ^ (Ly 

n 

(19) 

as the generators of an SU(2) isospin symmetry (orig- 
inally dubbed "axial charge" in Ref. l53l ). Since 
[-Himp-bathi I±] 0, the charge-coupled BFA model does 
not exhibit full isospin symmetry. However, this sym- 
metry turns out to be recovered in the asymptotic low- 
energy behavior at certain renormalization-group fixed 
points. 

As described in Ref. the computational effort re- 
quired for the NRG solution of a problem can be greatly 
reduced by taking advantage of these conserved quantum 
numbers. In particular, it is possible to obtain all phys- 
ical quantities of interest while working with a reduced 
basis of simultaneous eigenstates of S* 2 , S z , and Q with 
eigenvalues satisfying S z = S. With one exception noted 
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in Sec. IIV G[ any N s value specified below represents the 
number of retained (S, Q) multiplets, corresponding to a 
considerably larger number of (S, S z , Q) states. 

Even when advantage is taken of all conserved quan- 
tum numbers, NRG treatment of the charge-coupled 
BFA model remains much more demanding than that of 
the Anderson model [Eq. (fTTj) ] or the Anderson-Holstein 
model [Eq. (TlU|) ]. Being nondispersive, the bosons in the 
last model enter only the atomic- limit Hamiltonian Hq, 
allowing solution via the standard NRG iteration proce- 
dure. For Bose-Fermi models such as -ffccBFA, the need 
to extend a bosonic chain as well as a fermionic one at 
every even- numbered iteration N > 0, expands the ba- 
sis of Hm from AN S states to 4(Ab + 1)N S states, and 
multiplies the CPU time by a factor ~ (Nf, + l) 3 . Since 
we typically use iV b = 8 or 12 in our calculations, the 
increase in computational effort is considerable. 

The choice of value for the NRG discretization pa- 
rameter A involves trade-offs between discretization er- 
ror (minimized by taking A to be not much greater 
than 1) and truncation error (reduced by working with 
A ^> 1). Experience from other problem o 38 i 39 ' 52 indicates 
that critical exponents can be determined very accurately 
using quite a large A. Most of the results presented in the 
remaining sections of the paper were obtained for A = 9, 
with A — 3 being employed in the calculation of the im- 
purity spectral function. For convenience in displaying 
these results, we set ft = D = 1 and omit all factors of 
Po and Kq. 



III. PRELIMINARY ANALYSIS 

We begin by examining the special cases in which the 
impurity level is decoupled either from the conduction 
band or from the bosonic bath. Understanding these 
cases allows us to establish some expectations for the 
behavior of the full model described by Eq. ([T|). 



erators 



+ — == — (n d -l), 



>N q w q 



in the trivially solvable form 



ft' 

imp 



7 , W q a n li ,q ffi n d ,qi 



where 



ft' 

imp 



5 d (h d -l) + ^f(n d -l) 2 . 



(21) 



(22) 



(23) 



The bosons act on the impurity to reduce the Coulomb 
interaction from its bare value U to an effective value 



r,,, = c- — V^ 1 r - 



N„ 



B(u) 



dio. (24) 



For the bath spectral density in Eq. © with —1 < 
s < 0, one finds that for any nonzero e-b coupling A, 
Ucs — — °o and the singly occupied impurity states drop 
out of the problem. For the remainder of this section, 
however, we will instead focus on bath exponents s > 0, 
for which Eqs. © and ([Ml) give 



U, 



eff 



(25) 



For weak e-b couplings, U c g is positive and the ground 
state of i?zH lies in the sector n d = 1 where the impurity 
has a spin z component However, U e g is driven 

negative for sufficiently large A, placing the ground state 
in the sector n d — or n d — 2 where the impurity is 
spinless but has a charge (relative to half filling) of — 1 
or +1. 

Figure [1] illustrates this renormalization of the 
Coulomb interaction for the symmetric model (6 d = 0), 
in which the n d = and n d = 2 states always have the 
same energy. In this case, all four impurity states become 
degenerate at a crossover e-b coupling 



A. Zero hybridization 



K XcO = y/nsU /2Q,. 



(26) 



If one sets T = in Eq. (fT]), then the conduction band 
completely decouples from the remaining degrees of free- 
dom and can be dropped from the model, leaving the 
zero-hybridization model 



Hzh = S d {h d -l) + —(n d - l) 2 + ^w q a q a q 

q 

+ -^=(r ld -l)^A q (a q + a t _ q ). (20) 



The Fock space separates into sectors of fixed impurity 
occupancy {n d — 0, 1, or 2), within each of which the 
Hamiltonian can be recast, using displaced-oscillator op- 



The impurity contributions to physical properties at this 
special point, which is characterized by effective param- 
eters r = U = e d = 0, are identical to those at the 
free-orbital fixed point^ of the Anderson model. 

For the general case of an asymmetric impurity, the 
sectors n d — and 2 have a ground-state energy differ- 
ence Eo(n d = 2) — Eo(n d = 0) = 2S d for any value of A. 
The overall ground state of Eq. ([2H)l is a doublet (n d = 1 , 
S = ±i) for small e-b couplings, crossing over to a sin- 
glet (n d = for 8 d > 0, or n d = 2 for 5 d < 0) for large 
A. At AToAco = y / irs(U/2 — \5 d \)/£l, a point of three- fold 
ground-state degeneracy, the impurity contributions to 
low-temperature (T <ti \5 d \) physical properties are iden- 
tical to those at the valence-fluctuation fixed poinb 2 ^ of 
the Anderson model. 
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U_ 

2 



» % n d =0,2 



.2 

* ^ Aco 



* 1 



4- ft, 



A = 0. For s < 1, by contrast, the restriction &J„6 m < iV& 
leads, for A > A c o and large iteration numbers, to an ar- 
tificially truncated spectrum that cannot reliably access 
the low-energy physical properties. Nonetheless, obser- 
vation of this "localized" bosonic spectrum serves as a 
useful indicator, both in the zero-hybridization limit and 
in the full charge- coupled BFA model, that the effective 
e-b coupling remains nonzero. 

Another interpretation of Eq. (|27[) is that at the en- 
ergy scale E = Vth.~ k characteristic of interval k, the e-b 
coupling takes an effective value \(E) governed by the 
renormalization-group equation 



FIG. 1: Symmetric, zero-hybridization model defined by Hzh 
in Eq. (|2Up with 8d = 0: evolution with e-b coupling A 2 of the 
lowest eigenenergy in the spin sector (rid = 1> solid line) and 
in the charge sector (rid = 0, 2, dashed line). A level crossing 
occurs at A = A c o specified in Eq. (l26jl . 
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Using the NRG with only a bosonic chain [Eq. ([T3"f ] 
coupled to the impurity site, we have confirmed the ex- 
istence for Sd — of a simple level crossing from a spin- 
doublet ground state for A < A c0 to a charge-doublet 
ground state for A > A c o- In the former regime, the 
bosons couple only to the high-energy (rid — 0, 2) im- 
purity states, so the low-lying spectrum is that of free 
bosons obtained by diagonalizing given in Eq. 

(fT5)l . Here, NRG truncation plays a negligible role pro- 
vided that one works with Nb > 8 (say). 

For A > A c o, the low-lying bosonic excitations should, 
in principle, correspond to noninteracting displaced os- 
cillators having precisely the same spectrum as the orig- 
inal bath. However, the occupation number a q a q in 
the ground state of Eq. (|22[) obeys a Poisson distri- 



bution with mean \^/(N q u)q). Thus, the total num- 
ber of bosons corresponding to operators a q satisfying 
fiA~( fc+1 ) <w q < QA~ k takes a mean value 



(rik)o 



QA" 



B(u) 



OA-< fc + 1 > 

r (^oA) 2 



= < 



In A 



for s = I 



(A A) 2 (A 1 '" 1) 

7T (1-S) 



A (l-o)fc otherf 



(27) 



The bath states in the fcth interval are represented by 
NRG chain states < m < k, with the greatest weight 
being borne by state m — k. Thus, a faithful repre- 
sentation of the displaced-oscillator spectrum requires 
inclusion of states having b\ n b m up to several times 
(^m)o; based on experience with the Anderson-Holstein 
model^ one expects Nb > 4(ri m )o to suffice. Given that 
{n m )o oc A^ 1_s ^ m , it is feasible to meet this condition as 
m — > co so long as the bath exponent satisfies s > 1. 
Indeed, for Ohmic and super-Ohmic bath exponents, the 
NRG spectrum for A not too much greater than A c o is 
found to be numerically indistinguishable from that for 



d\n(n/E) 



(28) 



which implies that the e-b coupling is irrelevant for s > 
1, marginal for s = 1, and relevant for s < 1. While 
the NRG method is capable of faithfully reproducing the 
physics of -ffccBFA for arbitrary renormalizations of e^, 
U, and r, its validity is restricted to the region 



<- nN B 1-s a->i 
( k qA) ^ 



4 A 1 -* - 1 



4 In A 



(29) 



For A = 9 and Nb = 8, as used in most of our calcula- 
tions, the upper limit on the "safe" range of AToA varies 
from 1.7 for s = 1 to 0.9 for s — 0. 

We now focus on the value of the crossover e-b cou- 
pling A c o determined using the NRG approach. Figure 
[5] shows for five different bosonic bath exponents s that 
AoA c o has an almost linear dependence on the NRG dis- 
cretization A in the range 1.6 < A < 4. We believe that 
the rise in AoA c o with A reflects a reduction in the effec- 
tive value of Ko arising from the NRG discretization. It 
is knowr*2£ that in NRG calculations for fermionic prob- 
lems, the conduction-band density of states at the Fermi 
energy takes an effective value 



where 



p(0) = Po = p /A A , 



_ In A 1 + A" 1 

~ ~T~ i - Ta 31 



(30) 



(31) 



The general trend of the data in Fig. [5] is consistent with 
there being an analogous reduction of the bosonic bath 
spectral density that requires the replacement of K by 



K = K /A 



A,s 



(32) 



when extrapolating NRG results to the continuum limit 
A = 1. However, we have not obtained a closed- form 
expression for Aa, s - 

Table Q] lists values A c o(A — > 1) extrapolated from the 
data plotted in Fig. [21 For s > 0.4, these values are in 
good agreement with Eq. (|26j) . For s = 0.2, however, 
the extrapolated value of A c o lies significantly above the 
exact value, indicating that for given A the NRG un- 
derestimates the bosonic renormalization of U. This is 
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FIG. 2: Dependence of the level-crossing coupling A c o on the 
discretization A for the NRG solution of Hzh [Eq. with 
U = 0.1, 8 d = 0, N s = 200, iV& = 16, and five different values 
of the bath exponent s. Dashed lines show linear fits to the 
data. 



TABLE I: Crossover coupling A c0 for H Z h [Eq. J20j)] with 
U — 0.1, Sd = 0, and five different values of the bath ex- 
ponent s: Comparison between A c o(exact) given by Eq. (|26[) 
and A c o(A — * 1), the extrapolation to the continuum limit 
of numerical values obtained for N s == 200, Nb = 16, and 
1.6 < A < 4. Parentheses surround the estimated nonsystem- 
atic error in the last digit. 



s 


0.2 


0.4 


0.6 


0.8 


1.0 


A c o (exact) 


0.177 


0.251 


0.307 


0.355 


0.396 


A c0 (A^l) 


0.188(4) 


0.250(2) 


0.307(2) 


0.355(2) 


0.397(3) 



most likely another consequence of truncating the basis 
on each site of the bosonic tight-binding chain. 

In analyzing our NRG results for the full charge- 
coupled BFA model, we attempt to compensate for the 
effects of discretization and truncation by replacing Eq. 
(123) by 



77 m 



NRG 



U 



1 



(A/A c0 ) 2 



(33) 



Here, A c o is not the theoretical value predicted in Eq. 
(|26p . but rather is obtained from runs carried out for 
r = but otherwise using the same model and NRG 
parameters as the data that are being interpreted. 



B. Zero electron-boson coupling 

For A = 0, the bosonic bath decouples from the elec- 
tronic degrees of freedom, which are then described by 
the pure Anderson model. In this section, we briefly 
review aspects of the Anderson model that will prove 
important in interpreting results for the charge-coupled 
BFA model. For further details concerning the Anderson 
model, see Refs. [l| and 120. 



For any T > 0, and for any U and Sd = £d + U/2 
(whether positive, negative, or zero), the stable low- 
temperature regime of the Anderson model lies on a line 
of strong- coupling fixed points corresponding to T = oo. 
At any of these fixed points, the system is locked into the 
ground state of the atomic Hamiltonian Ho, and there are 
no residual degrees of freedom on the impurity site or on 
site n = of the fermionic chain; the NRG excitation 
spectrum is that of the Hamiltonian 20 



(34) 



The coefficients r„ are identical to those entering H^^9 



[Eq. (H])], except that here n = 0. Note that in Eq. 
the sum over n begins at 1 rather than 0. 

As shown in Ref. [2(| the strong-coupling fixed points 
of the Anderson model are equivalent — apart from a shift 
of 1 in the ground-state charge Q defined in Eq. (fT5)l — 
to the line of frozen-impurity fixed points corresponding 
to €d — oo, r = U = 0, with NRG excitation spectra 
described by 



NRG 
FI 



(35) 



The mapping between alternative specifications of the 
same fixed-point spectrum is^°- 



TjOoVo = ~(np Vi)~ 



(36) 



where po [see Eq. (|3U)) ] is the effective conduction-band 
density of states. 

The fixed-point potential scattering is related to the 
ground-state impurity charge via the Friedel sum rule, 

2 2 
(fid— l)o = — arccot(7rp Vo) = — arctanf— npoVi). (37) 

For |<Jd|,r one finds that 



(fid - l)o = - 



ttA a U 2 ' 



(38) 



where A\ is defined in Eq. (|3l|) . 

Even though the stable fixed point of the Anderson 
model for any T > is one of the strong-coupling fixed 
points described above, the route by which such a fixed 
point is reached can vary widely, depending on the rela- 
tive values of U, Sd, and T. For our immediate purposes, 
it suffices to focus on the symmetric (Sd — 0) model, for 
which there is a single strong-coupling fixed point corre- 
sponding to Vq = ±oo or V\ = 0. If the on-site Coulomb 
repulsion is strong enough that the system enters the 
local-moment regime (T,T <C U), then it is possible to 
perform a Schrieffer- Wolff transformation^ that restricts 
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the system to the sector n<j = 1 and reduces the Anderson 
model to the Kondo model described by the Hamiltonian 

H K = flband+J^- ("dt - n d[ ) E (4t c k'T - c ki c k'i 
k k,k' 



performing a Lang-Firsov^ transformation TJccbfa 



^CCBl \ 



[/ = exp 



[/ 1 i?ccBFA^ with 
A 



fa - 1) E ■ 



= (a„ - aj.) 



(43) 



2jy fc E] (^|^i c k| c k'T H.c.^ , (39) The transformation eliminates Hi mp -bath> leaving 

k,k' 



where 



sr 

PoJ± = —fj- 

TTU 



(40) 



The stable fixed point is approached below an exponen- 
tially small Kondo temperature 7k when the spin-flip 
processes associated with the J± term in 77k cause the 
effective values of poJ z and PqJ± to renormalize to strong 
coupling, resulting in many-body screening of the impu- 
rity spin. 

Motivated by the discussion in Sec. IIII A[ we also con- 
sider the case of strong on-site Coulomb attraction. In 
the local-charge regime (T, T <C —U), a canonical trans- 
formation similar to the Schrieffer- Wolff transformation 
restricts the system to the sectors n d = and rid = 2, and 
maps the Anderson model onto a charge Kondo model 
described by the Hamiltonian 

H CK = ffband + {n d - 1) E (4 T c k , T + 4i«V; 



+ iHr E { d \ d Wi c ^ + H - c -) . ( 41 ) 

k k,k' 



k,k' 



where 



PoW d = p W p = 



2T 

W\' 



(42) 



In this case, the stable fixed point is approached be- 
low an exponentially small (charge) Kondo temperature 
Tk when the charge-transfer processes associated with 
the W p term in Hck cause the effective values of poWd 
and PqW p to renormalize to strong coupling, resulting 
in many-body screening of the impurity isospin degree 
of freedom [associated with the d-operator terms in Eqs. 

Between the opposite extremes of large positive U and 
large negative U is a mixed- valence regime T, \ U\ <C T 
in which interactions play only a minor role. Here, the 
stable fixed point is approached belo w a tempe rature of 
order F when the effective value of y/T/ (2itD) scales to 
strong coupling, signaling strong mixing of the impurity 
levels with the single-particle states of the conduction 
band. 



C. Expectations for the full model 

Insight into the behavior of the full charge-coupled 
BFA model described by Eqs. ([I])-© can be gained by 



^CCBFA — -^imp + -^band 



-ffbath + -^imp-band' 



(44) 



where H[ mp is as defined in Eqs. (|2"3"|) and (|2"4")l . and 



imp-band 



Es^exp 



E 



A q (e 



V£ exp 



E 



A q (a q - a q ) 



'N q w q 



(45) 



In addition to renormalizing the impurity interaction 
from U to U c b entering 77;' mp , the e-b coupling causes 
every hybridization event to be accompanied by the cre- 
ation and annihilation of arbitrarily large numbers of 
bosons. 

In the case of the Anderson- Holstein model [Eq. (|T0|) ]. 
various limiting behaviors are understood^ In the in- 
stantaneous limit ujq 3> r, the bosons adjust rapidly to 
any change in the impurity occupancy; for Xq/ojo <C U <C 
wo, the physics is essentially that of the Anderson model 
with U — > U e s, while for Xq/uj 3> D, U,T, there is also 
a reduction from T to Texp[— (Ao/wo) 2 ] in the rate of 
scattering between the n d = and nd — 2 sectors, re- 
flecting the reduced overlap between the ground states 
in these two sectors. In the adiabatic limit ojq <C T, the 
phonons are unable to adjust on the typical time scale 
of hybridization events, and neither U nor T undergoes 
significant renormalization. 

Similar analysis for the charge- coupled BFA model is 
complicated by the presence of a continuum of bosonic 
mode energies w, only some of which fall in the instanta- 
neous or adiabatic limits. Nonetheless, we can use results 
for the cases T = (Sec. [TO} and A = (Sec.[mB|, as 
well as those for the Anderson-Holstein model, to identify 
likely behaviors of the full model. Specifically, we focus 
here on the evolution with decreasing temperature of the 
effective Hamiltonian describing the essential physics of 
the symmetric (e^ = —U/2) model at the current tem- 
perature. This effective Hamiltonian is obtained under 
the assumption that real excitations of energy above the 
ground state E > rfT — where 77 is a number around 5, 
say — make a negligible contribution to the observable 
properties, and thus can be integrated from the problem. 

Based on the preceding discussion, one expects that 
at high temperatures T ^> T, the physics of the charge- 
coupled BFA model will be very similar to that of the 
Anderson model with U replaced by U(rjT), where 



U(E) = U - - 



B(co) 



diu. 



(46) 
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Note that (7(0) is identical to U c g defined in Eq. (|24f . 
For the bath spectral density in Eq. ([5]) with s > 0, 

U(E) = U- 2(X ° A)2 [i _ (E/QY]Q. (47) 

TTS 

When analyzing NRG data, we instead use 

U miG (E) = U{\ - (A/A c0 ) 2 [1 - (E/n) s ] }, (48) 

where A c o is the empirically determined value discussed 
in connection with Eq. (|33l) . 

If, upon decreasing the temperature to some value 
Tlm, the system comes to satisfy U^T^m) = 
t] max(TLM j r) , then it should enter a local-moment 
regime described by the effective Hamiltonian H^m = 
Hk + -ffbath with the exchange couplings in Hk [Eq. (f3"9")) ] 
determined by Eq. (jlO")) with U — > U{riT LM ), similar to 
what is found in the Anderson-Holstein model. 8 Since 
they couple only to the high-energy sectors n& = and 
iid = 2 that are projected out during the Schrieffer- Wolff 
transformation, the bosons should play little further role 
in determining the low-energy impurity physics. The out- 
come should be a conventional Kondo effect where the 
e-b coupling contributes only to a renormalization of the 
Kondo scale Tk. 

If, instead, at some T = Tlc the system satisfies 
C/(?/Tlc) = — i] max(XLCj T), then it should enter a local- 
charge regime described by the effective Hamiltonian 

HhC — HcK + #bath + -^imp-bath- (49) 

Based on the behavior of the Anderson-Holstein modelj 8 
one expects Wd in Hqk [Eq. (|4Tj) ] to be determined by 
Eq. flUD with U -> U(r)T LC ), but with W p exponentially 
depressed due to the aforementioned reduction in the 
overlap between the ground states of the = and 
rid = 2 sectors. The bosons couple to the low-energy sec- 
tor of the impurity Fock space, so they have the potential 
to significantly affect the renormalization of Wd and W p 
upon further reduction in the temperature. In particu- 
lar, the A term in H^c which favors localization of the 
impurity in a state of well-defined rid = or 2, directly 
competes with the W p double-charge transfer term that is 
responsible for the charge Kondo effect of the negative-?/ 
Anderson model. This nontrivial competition gives rise 
to the possibility of a QPT between qualitatively distinct 
ground states of the charge-coupled BFA model. 

Between these extremes, the system can enter a mixed- 
valence regime of small effective on-site interaction. In 
this regime, one must retain all the impurity degrees 
of freedom of the charged-coupled BFA model. The 
impurity-band hybridization competes with the e-b cou- 
pling for control of the impurity, again suggesting the 
possibility of a QPT. 

Each of the regimes discussed above features competi- 
tion between band-mediated tunneling within the man- 
ifold of impurity states and the localizing effect of the 
bosonic bath. Although the tunneling is dominated by a 



different process in the three regimes, it always drives the 
system towards a nondegenerate impurity ground state, 
whereas the e-b coupling favors a doubly-degenerate 
(rid = 0, 2) impurity ground state. In order to provide 
a unified picture of the three regimes (and the regions of 
the parameter space that lie in between them), we will 
find it useful to interpret our NRG result in terms of an 
overall tunneling rate A, which has a bare value 

A ~ +2TD/tt+ (50) 

Here, W p is assumed to be negligibly small in the local- 
moment regime, and J± to be similarly negligible in the 
local-charge regime. If A renormalizes to large values 
while the e-b coupling A scales to weak coupling, then 
one expects to recover the strong-coupling physics of the 
Anderson model. If, on the other hand, A becomes strong 
while A becomes weak, the system should enter a low- 
energy regime in which the bath governs the asymptotic 
low-energy, long-time impurity dynamics. Whether or 
not each of these scenarios is realized in practice, and 
whether or not there are any other possible ground states 
of the model, can be determined only by more detailed 
study. These questions are answered by the NRG results 
reported in the sections that follow. 

IV. RESULTS: SYMMETRIC MODEL WITH 
SUB-OHMIC DISSIPATION 

This section presents results for Hamiltonian (p} with 
U = —led > and with sub-Ohmic dissipation char- 
acterized by an exponent < s < 1 . Figure [3J shows 
a schematic phase diagram on the A-r plane at fixed U. 
There are two stable phases: the localized phase, in which 
the impurity dynamics are controlled by the coupling to 
the bosonic bath and the system has a pair of ground 
states related to one another by a particle-hole transfor- 
mation; and the Kondo phase, in which there is a nonde- 
generate ground state. These phases are separated by a 
continuous QPT that takes place on the phase boundary 
(solid line in Fig. [3]), which we parametrize as A = A c (r). 
Within the Kondo phase, the nature of the correlations 
evolves continuously with increasing A (at fixed T) from 
a pure spin-Kondo effect for A = to a predominantly 
charge-Kondo effect beyond a crossover (dashed line in 
Fig. [3]) associated with the change in sign of U s defined 
in Eq. (24]). 

As s decreases, and the e-b coupling becomes increas- 
ingly relevant — in a renormalization-group sense [see Eq. 
(|28|] — the phase boundary moves to the left as the lo- 
calized phase grows at the expense of the Kondo phase, 
which disappears entirely for s < 0. As will be seen in 
Sec. [V] the phase diagram of the Ohmic (s — 1) problem 
has the same topology as Fig. [3J even though (as de- 
scribed in Sec. [Vj the nature of the QPT is qualitatively 
different than for < s < 1. For s > 1, the e-b coupling 
is irrelevant, and the system is in the Kondo phase for 
all T > 0. 
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FIG. 3: Schematic phase diagram of the symmetric charge- 
coupled BFA model for bath exponents < s < 1. The solid 
curve marks the boundary between the Kondo phase, in which 
the impurity degrees of freedom are screened by conduction 
electrons, and the localized phase, in which the impurity dy- 
namics are controlled by the coupling to the dissipative bath. 
The dashed vertical line represents a crossover from a regime 
in which Kondo screening takes place primarily in the spin 
sector to a regime in which a charge-Kondo effect is predom- 
inant. 



The remainder of this section presents the evidence 
for the previous statements. We first discuss the 
renormalization-group flows and fixed points. We then 
turn to the behavior in the vicinity of the phase bound- 
ary, focusing in particular on the critical response of the 
impurity charge to a local electric potential. Following 
that, we present results for the impurity spectral func- 
tion, and show that the low-energy scale extracted from 
this spectral function supports the qualitative picture 
laid out in the paragraphs above and summarized in Fig. 

m 



A. NRG flows and fixed points 

Figure 0] plots the schematic renormalization-group 
flows of the couplings A entering Eq. (| 1 5|) and A de- 
fined in Eq. (|50[) for a symmetric impurity (U — — 2e^) 
coupled to bath described by an exponent < s < 1. 
These flows are deduced from the evolution of the many- 
body spectrum with increasing iteration number N, i.e., 
with reduction in the effective band and bath cutoffs 
D = Q ~ DA~ N / 2 . A separatrix (dashed line) forms 
the boundary between the basins of attraction of a pair 
of stable fixed points, regions that correspond to the two 
phases shown in Fig. [3l Figure [4] also shows three un- 
stable fixed points. In contrast to the situation at other 
points on the flow diagram, each of the fixed points ex- 
hibits a many-body spectrum that can be interpreted as 
the direct product of a set of bosonic excitations and a 
set of fermionic excitations. 

The Kondo fixed point corresponds in the 
renormalization-group language of Fig. 2] to effec- 
tive couplings A = and A = oo. The many-body 




FIG. 4: Schematic renormalization-group flows on the A-A 
plane for the symmetric charge-coupled BFA model with a 
bath exponent < s < 1. Trajectories with arrows represent 
the flow of the couplings A entering Eq. (|15[) and A defined 
in Eq. (150[l under decrease of the high-energy cutoffs on the 
conduction band and the bosonic bath. Between the basins of 
attraction of the Kondo fixed point (K) and the localized fixed 
point (L) lies a separatrix, along which the flow is away from 
the free-orbital fixed point (FO) located at A = A c q, A = 
and toward the critical fixed point (C). For A = only, there 
is flow from FO towards the local-moment fixed point (LM) 
at A = 0. 

spectrum decomposes into the direct product of (i) the 
excitations of a free bosonic chain described by Eq. Q13]) 
alone, and (ii) the strong-coupling excitations of the 
Kondo (or symmetric Anderson) model, corresponding 
to free electrons with a Fermi- level phase shift of tt/2. 
This spectrum, which exhibits SU(2) symmetry both 
in the spin and charge (isospin) sectors, is identical 
to that found throughout the Kondo phase of the 
particle-hole-symmetric Ising BFK Hamiltonia n 38 i 39 (a 
model in which the bosons couple to the impurity's spin 
rather than its charge). 

The schematic RG flow diagram in Fig. [4] shows a lo- 
calized fixed point corresponding to A = oo and A = 0. 
However, this is really a line of fixed points described by 
Hlc [Eq- (HU)] with effective couplings A = oo, W p = 0, 
and < Wd < oo. Since W p = 0, the impurity occu- 
pancy takes a fixed value nd — or 2. (It is important to 
distinguish rid, used to characterize the fixed-point exci- 
tations, from the physical expectation value of fid- The 
latter quantity is discussed in Sec. IIVFTT1 ) 

Each fixed point along the localized line has an excita- 
tion spectrum that decomposes into the direct product of 
(i) bosonic excitations identical to those at the localized 
fixed point of the spin-boson models with the same bath 
exponent s, and (ii) fermionic excitations described by a 
Hamiltonian 

HIT = « + W d {n d - 1) (T flf 0a , (51) 

which is just the discretized version of Hck [Eq. (HU] 
with W p = and the operator fid replaced by the param- 
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FIG. 5: (Color online) Fixed-point coupling Wd entering Eq. 
(|51l) vs e-b coupling A — A c in the localized phase near the 
phase boundary at A = A c . Results are shown for U = — 2ed = 
0.1, A = 9, N s = 500, N b = 8, four different values of the bath 
exponent s, and T = 0.5, 1.0, 10, and 50 for s = 0.2, 0.4, 0.6, 
and 0.8, respectively (Ref. |56j). The power-law divergence of 
Wd as A — > A^ [Eq. (|52[)] is reflected in the linear behaviors 
of data on a logarithmic scale. The numerical values of the 
exponent /3 obtained here are identical (to within small errors) 
to those listed in Table ITlTI 



eter rid- The low- lying many-body eigenstates of H^t 
appear in degenerate pairs, one member of each pair cor- 
responding to rid = and the other to rid = 2. The 
fixed-point coupling Wd increases monotonically as the 
bare e-b coupling A decreases from infinity, and diverges 
on approach to the phase boundary. As illustrated in Fig. 
[H this divergence can be fitted to the power-law form 



W d oc (A - \ c )- p 



for A -> Aj 



(52) 



For reasons that will be explained in Sec. IIVE 11 the nu- 
merical value of /3 coincides, to within a small error, with 
that of the order-parameter exponent [3 defined in Eq. 

CD- 

The free-orbital fixed point (A = A c o, A = 0) is unsta- 
ble with respect to a bare T ^ or any deviation of A 
from A c o = limr^oA c (r). The local-moment fixed point 
(A = A = 0), at which the impurity has a spin-i degree 
of freedom decoupled from the band and from the bath, 
is reached only for bare couplings T = (hence, A = 0) 
and A < A c o. 

Of greatest interest is the unstable critical fixed point 
that is reached for any bare couplings lying on the bound- 
ary A = A c (r) between the Kondo and localized phases. 
At this fixed point, the low- lying spectrum can be con- 
structed as the direct product of (i) the critical spectrum 
of the spin-boson model with the same bath exponent s, 
and (ii) the strong-coupling spectrum of the Kondo (or 
symmetric Anderson) model. This spectrum, which ex- 
hibits full SU(2) symmetry in both the spin and isospin 
sectors, is identical to that at the critical point of the 
Ising-anisotropic Bose-Fermi Kondo model,— and is il- 
lustrated in Fig. 3(c) of Ref.[H. 



The decomposition of the critical spectrum can be un- 
derstood by considering the flow of couplings entering 
the local-charge Hamiltonian H^c defined in Eq. (|49[) . 
The fixed-point value of the density-density coupling is 
Wd = oo in the charge-Kondo regime of the Kondo phase 
and diverges according to Eq. ([52")) in the localized phase. 
It is therefore reasonable to assume that in the vicinity of 
the phase boundary, Wd rapidly renormalizes to strong 
coupling, locking the impurity site and site n — of 
the fermionic chain into one of just two states, which 
we can write in a pseudospin notation as Iff) = djd||0) 

and |4) = /o|/oj0), where |0) is the no-particle vacuum. 
Hopping of electrons on or off site n = is forbidden, 
so the discretized form of TJlc reduces to an effective 
Hamiltonian 



-Hi 



NRG 
LC 



(W d = oo) 



zYNRG 



(0) 



o-NRG 
-"SBM • 



(53) 



Here, HgQ G (0) [Eq. (f34|) ] acts only on fermionic chain 
sites n > 1, and yields the Kondo/ Anderson strong- 
coupling excitation spectrum, while 



NRG 
SBM 



= H 



NRG 



2W r 



bath Tf"f, 

VLK \ 



(b Q + bl) (54) 



acts on the remaining degrees of freedom in the problem 
in a subspace of states all carrying quantum numbers 
S = S z = Q = 0. H ggjjj is precisely the discretized 
form of the spin-boson Hamiltonian with tunneling rate 
A = AW P and dissipation strength a = 2(Ko\) 2 /n. 
These two couplings compete with one another, with 
three possible outcomes: (1) A can scale to infinity 
and a to zero, resulting in flow to the delocalized fixed 
point (the Kondo fixed point of the charge-coupled BFA 
model); (2) a can scale to infinity and A to zero, yield- 
ing flow to the localized fixed point; or (3) both couplings 
can renormalizc to finite values A = Ac, a — ac at the 
critical point. This picture implies that the universal crit- 
ical behavior of the charge-coupled BFA model should be 
identical to that of the spin-boson model, the conduction- 
band electrons serving only to dress the rid = 0, 2 impu- 
rity levels and to renormalize the impurity tunneling rate 
and the dissipation strength. 

Given that the NRG approach necessarily involves 
Fock-space truncation, it is instructive to examine the 
dependence of the fixed-point spectra on the parame- 
ters N s and Nb denoting, respectively, the number of 
states retained from one NRG iteration to the next and 
the maximum number of bosons allowed per site of the 
bosonic chain. Figure [6] shows, for representative bath 
exponents s = 0.2 and s — 0.8, that the energy of the 
lowest bosonic excitation at A = A c converges rapidly 
with increasing N s and JV&. This behavior suggests that 
for A = 9, at least, N s = 500 and Nj, = 8 are sufficient 
for studying the physics at the critical point. 

By contrast, the lowest bosonic excitation energy for 
A = 1.1 A c , plotted in Fig. [71 converges only slowly with 
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FIG. 6: Dependence of the energy of the first bosonic excita- 
tion at the critical point (A = A c ) on the NRG truncation pa- 
rameters iVf, and N a . Results are shown for U = — 2ej = 0.1, 
r = 0.01, A = 9, and bath exponents s = 0.2 and s = 0.8. In 
the left panels, N s — 500, while in the right panels Nt = 8. 



respect to N^. This points to the failure of the truncated 
bosonic basis deep inside the localized phase of the sub- 
Ohmic model, where the mean boson number per site is 
expected to diverge according to Eq. f2"7]). This inter- 
pretation is confirmed by calculation of the expectation 
value of the total boson number, 



B 



N 



M(N) 



(55) 



where M(N) denotes the highest labeled bosonic site 
present at iteration N. Our results for (-B20) vs Nh (not 
shown) are very similar to those in Fig. 8 of Ref. [39|, with 
convergence by Nf, = 8 at the critical point, but no evi- 
dence of such convergence for an e-b coupling 10% over 
the critical value. 

Recently, Bulla et al. applied a "star" reformulation of 
the NRG to the spin-boson model^ While this approach 
provides a good description of the localized fixed point, it 
does not correctly capture the physics of the delocalized 
phase (corresponding to the Kondo phase of the present 
model) or of the critical point that separates the two 
stable phases. For this reason, we prefer to work with 
the "chain" formulation summarized in Sec. Hfl 



B. Critical coupling 

Figure [5] plots the critical e-b coupling X C (T) for fixed 
U = —2ed and four different values of the bath exponent 
s. As expected, with increasing T, the critical coupling 
increases smoothly from A c (r = 0) = A c o, reflecting the 
fact that entry to the localized phase requires an e-b cou- 
pling sufficiently large not only to drive J7 e ff negative, but 
also to overcome the reduction in the electronic energy 
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FIG. 7: Dependence of the energy of the first bosonic excita- 
tion in the localized phase (A = 1.1A C ) on the NRG truncation 
parameters Nb and N s . All other parameters are as in Fig. [6] 



0.10 




FIG. 8: (Color online) Critical coupling A c vs hybridization 
width r for U = -2e d = 0.1, A = 9, N s = 500, N b = 8, and 
the bath exponents s listed in the legend. 



that derives from the hybridization. We believe that the 
vertical slope of the s = 0.2 phase boundary as it ap- 
proaches the horizontal axis in Fig. [8] is an artifact stem- 
ming from the same source as the NRG overestimate of 
A c o for the same bath exponent. (See the discussion of 
Fig. [Din Sec. UTS]) 

In the subsections that follow, we show that the critical 
properties of the charge-coupled BFA model map, under 
interchange of spin and charge degrees of freedom, onto 
those of the spin-coupled BFA model studied (along with 
the corresponding Ising BFK model) in Ref. [39, The 
spin-coupled model is described by Eqs. U)-© and (fl2|) - 
with Eqs. © and (T5]) replaced by 



imp-bath 



-.{n d] -n di )^g n 



(56) 
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FIG. 9: Variation with bath exponent s of the critical cou- 
plings A c and A c o in the charge-coupled BFA model (this 
work) and g c /2 in the spin-coupled BFA model (Ref. |39T ). 
Results are shown for U = ~2e d = 0.1, T = 0.01, A = 9, 
N s = 500, and N b = 8. 
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FIG. 10: (Color online) Crossover scale T* vs A c — A on the 
Kondo side of the critical point for four different values of the 
bath exponent s, with all other parameters as in Fig. [5] The 
slope of each line on this log-log plot gives the correlation- 
length exponent v(s) defined in Eq. (|58[) , 



and 



rVNRG 
imp-bath 



ClK g 



2vM^+l) 



(n d T ~ndi)( b o 



bl) 



(57) 



In light of the parallels between the universal critical be- 
havior of the two models, it is of interest to compare 
their critical couplings, making due allowance for the ad- 
ditional prefactor of \ that enters Eqs. (|56[) and ([571 



Figure [9] plots the s dependence of A c and g c /2 for 
fixed values of U — —2ed and T. For all < s < 1, 
A c is found to exceed g c /2- This fact can be understood 
by noting the contrasting role of the e-b coupling in the 
two models. In the spin-coupled BFA model, increasing 
g from zero immediately begins to localize the impurity 
in a state of fixed S z , and thereby to impede the spin- 
flip processes that are central to the Kondo effect. In the 
charge-coupled model, by contrast, increasing A from zero 
initially acts to decrease the effective Coulomb repulsion 
and hence to enhance charge fluctuations on the impurity 
site; only for A c > A c o do further increases in the e- 
b coupling serve to localize the impurity in a state of 
fixed charge, eventually leading to the suppression of the 
charge Kondo effect at A = A c . 



C. Crossover scale 

Under the renormalization-group flows sketched in Fig. 
[U the system passes, with decreasing energy cutoff or 
decreasing temperature, between the regions of influence 
of different renormalization-group fixed points. For bare 
parameters that place the system near the boundary be- 
tween the Kondo and localized phases, the free-orbital 
fixed point typically governs the behavior at tempera- 
tures much greater than the Kondo temperature 7k of 
the Anderson model obtained by setting A = in Eq. fl} . 
For temperatures between of order 7k and a crossover 



scale T*, the system exhibits quantum critical behav- 
ior controlled by thermal fluctuations about the unstable 
critical point. Finally, the physics in the regime T < 7 1 * 
is governed by one or other of the two stable fixed points: 
Kondo or localized. 

For fixed values of all other parameters, one expects 
T* to vanish as the e-b coupling approaches its critical 
value according to a power law: 



T* oc | A — \ c \ v for A — * A c 



(58) 



where v is the correlation-length exponent. 4 The 
crossover scale can be determined directly from the NRG 
solution via the condition T» cx A _Ar */ 2 , where is 
the number of the iteration at which the many-body en- 
ergy levels cross over to those of a stable fixed point. 
There is some arbitrariness as to what precisely consti- 
tutes crossover of the levels. Different criteria will pro- 
duce T*(A) values that differ from one another by a A- 
independent multiplicative factor. It is of little impor- 
tance what definition of iV» one uses, provided that it is 
applied consistently. 

Figure rTU] shows typical dependences of T* on A c — A 
in the Kondo phase. Equation ([58]) holds very well over 
several decades, as demonstrated by the linear behavior 
of the data on a log-log plot. We find that the numer- 
ical values of v(s), some of which are listed in Table HT1 
are identical (within small errors), to those of the spin- 
boson and Ising BFK models for the same bath exponent 
s. This supports the notion that the critical point of the 
charge-coupled BFA model belongs to the same univer- 
sality class as the critical points of the spin-boson and 
Ising BFK models. However, to confirm this equivalence, 
we must compare other critical exponents, as reported 
below. 
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TABLE II: Correlation-length critical exponent v vs bath ex- 
ponent s for the charge-coupled Bose-Fermi Anderson model 
(CC-BFA, this work) and for the Ising-anisotropic Bose-Fermi 
Kondo model (BFK, from Refs. l38l andl39T). Parentheses sur- 
round the estimated nonsystematic error in the last digit. 



s 


0.2 


0.4 


0.6 


0.8 


^(CC-BFA) 


4.99(3) 


2.52(2) 


1.97(4) 


2.12(6) 


v(BFK) 


4.99(5) 


2.50(1) 


1.98(3) 


2.11(2) 



D. Thermodynamic susceptibilities 

In this subsection, we consider the response of the 
charge-coupled BFA model to a global magnetic field 
H and to a global electric potential $. These external 
probes enter the Hamiltonian through an additional term 



H, 



HS Z + *Q, 



(59) 



where S z and Q are defined in Eqs. (| 16[) and (|T8|) , respec- 
tively. In particular, we focus on the static impurity spin 
susceptibility Xs,im P = —d 2 Fi mp /dH 2 and the static im- 
purity charge susceptibility Xc,imp = <9 2 F imp /<9$ 2 . Here, 
F lmp = A(F), where A(X) is the difference between (i) 
the value of the bulk property X when the impurity is 
present and (ii) the value of X when the impurity is re- 
moved from the system. It is straightforward to show 
that 



r Xs , imp = a(((£ 2 2 )) - «&)) 2 ; 

T Xc , imp = A(((Q 2 )> - ((Q)> 2 ) 
where, for any operator A, 

Triexp(-F/T) 



((A)) = 



Tr exp(-ff/T) 



(60) 
(61) 



(62) 



Note that with the above definitions, limT^oo Tx s ,imp = 
| but limj^oo T^cimp = \-, a factor of four differ- 
ence that must be taken into account when comparing 
the two susceptibilities. Since each Tximp is calculated 
as the difference of bulk quantities, its evaluation us- 
ing the NRG method is complicated by significant dis- 
cretization and truncation errors. In order to obtain 
reasonably well-converged results for Tximp, we retain 
N s = 2000 states after each NRG iteration. However, 
even this number is insufficient to allow reliable extrac- 
tion of Ximp = (Txim P )/T as T -> 0. 

Figure [ft] plots NRG results for T Xs ,im P {T) and 
jTxc,imp{T), calculated for bath exponent s = 0.8 and 
different values of the e-b coupling A. For A <C A c o (see 
Sec. IIII Aj) . both impurity susceptibilities behave very 
much as they do in the Anderson model: with decreasing 
temperature, Tx c , imp quickly falls toward zero, signal- 
ing quenching of charge fluctuations upon entry into the 
local-moment regime, whereas Tx s ,imp initially rises to- 
wards its local- moment value of 4, before dropping to 



zero for T <C T* on approach to the Kondo fixed point. 
With increasing A, the charge response grows and the 
spin response is suppressed. The two susceptibilities are 
approximately equivalent for A — A c o, where the effec- 
tive Coulomb interaction U c s — 0. For still stronger e-b 
couplings, Txs,imp plunges rapidly as the temperature 
is decreased, whereas Tx c , imp first rises on entry to the 
local-charge regime before dropping to satisfy 



Urn Txc,im P {T) = for A < A c 



(63) 



These trends are very similar to those exhibited-^ by the 
Anderson-Holstein model. In that model, however, the 
drop in Tx c , imp{T) takes placed for strong e-b couplings 
Ao 3> \/u)qU/2 around an effective Kondo temperature 
Tf ~ £>exp(-vrA(5/rwg). In the charge-coupled BFA 
model, by contrast, neither the spin susceptibility nor the 
charge susceptibility exhibits any obvious feature that 
correlates with the vanishing of T» as A — > X~ . This can 
be understood by noting that the impurity susceptibili- 
ties are determined purely by the fermionic part of the 
excitation spectrum, whose asymptotic low-energy form 
is the same at the critical fixed point (which governs the 
behavior in the quantum critical regime T* < T < Tr-) 
as at the Kondo fixed point (which controls the regime 
T < T 1 *). 

The behavior of the static impurity spin susceptibil- 
ity is qualitatively unchanged upon crossing from the 
Kondo phase to the localized phase. However, for A > 
A c , Txc,im P approaches at low temperatures a nonzero 
value that can be inferred from the effective Hamiltonian 
ijNRG [Eq. I|5ip]. Electrons near the Fermi level experi- 
ence an s-wave phase shift 



S(u = 0) 



So for rid = 

7T - S for n d = 2, 



(64) 



where rid labels the two disconnected sectors of H^^P, 
and 



Sq = arctan (npoWd), < So < tt/2, 



(65) 



with po being the effective conduction-band density of 
states defined in Eq. (|30|) . It is then straightforward to 
show that 



HmTxctopCT) = (l-2,5o/7r) 2 



(66) 



Equations (|52"]) . (|65[) . and (|66[) together imply that 

HmT X c,ta.p(T) « (X-Xc) 20 for A -> A+ (67) 

As this example illustrates, the thermodynamic sus- 
ceptibilities contain signatures of an evolution from a 
spin-Kondo effect to a charge-Kondo effect. Further- 
more, Eqs. (|6"3"|) and (f6"T)) suggest that Xc,im P may serve 
as the order-parameter susceptibility for the QPT. How- 
ever, neither susceptibility manifests the vanishing of the 
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FIG. 11: (Color online) Temperature dependence of the im- 
purity contribution to the static spin (left) and charge (right) 
susceptibilities for s = 0.8, U = ~2e d = 0.1, T = 0.01, 
A = 9, N s = 2000, N b = 8, and different values of the e-b 
coupling A. Dotted curves correspond to e-b couplings lying 
between the A values specified in the legend for the adjacent 
nondotted curves. For A = A c o — 0.396, the spin and charge 
susceptibilities are equivalent: Xa,im P {T) ~ jXc,im P (T). For 
A < A c o, the spin response is stronger, while for A > A c o, 
the charge response dominates. For A < A c ~ 0.5052181, 
liniT-vQ Txc,imp(T) = 0, whereas for A > A c , the limiting 
value is nonzero and obeys Eqs. (|65[l and (|66[l . 



crossover scale T* on approach to the transition from the 
Kondo side. Moreover, the conservation of Q prevents 
Xc,imp from acquiring an anomalous temperature depen- 
dence in the quantum-critical regime^ Thus, one is led 
to conclude that the response to a global electric poten- 
tial $ does not provide access to the critical fluctuations 
near the QPT. 



E. Local charge response 

Given the nature of the coupling in Hamiltonian ((T|) 
between the impurity and the bosonic bath, we expect to 
be able to probe the quantum critical point through the 
system's response to a local electric potential 4> that acts 
solely on the impurity charge, entering the Hamiltonian 
via an additional term 



-ff c ,ioc = (j>{n d - I). 



(68) 



A nonzero <fi is equivalent to a shift in Sd entering Eq. §§§ 
away from its bare value ed + U/2 = 0. 

In this subsection we show that for sub-Ohmic bath 
exponents < s < 1, (i) the response to a static 4> is 
described by critical exponents that satisfy hyperscaling 
relations characteristic of an interacting quantum critical 
point, (ii) numerical values of these critical exponents 
are identical to those of the spin-boson and Ising BFK 
models, and (iii) the dynamical response is consistent 



with the presence of ui/T scaling in the vicinity of the 
quantum critical point. 



1. Static local charge response 

The response to imposition of a static local potential <f> 
is measured by the thermodynamic average value of the 
impurity charge, 

Qloc = {{fid - 1», (69) 
and through the static local charge susceptibility 

T ^loc 



Xc,ioc(T; uj = 0) = 



= — lim 

0=o 0-0 4) 



(70) 



In NRG calculations of lim^o Qioc(</ ) ) and Xc.ioc, we use 
potentials in the range 10 -13 < \(j>\ < 10 -10 . 

As illustrated in Fig. 1121 the "spontaneous impurity 
charge" lim^o Qi oc (A, (f>; T — 0) indeed serves as an or- 
der parameter for the QPT between the Kondo and lo- 
calized phases. This quantity vanishes for all A < A c and 
is nonzero for A > A c , its onset being described by the 
power law 

limQ loc (A,0;T = O)cx(A-A c ) /3 for A -> A+. (71) 

In the localized phase, the presence of an infinitesimal 
local potential restricts the effective Hamiltonian (|5Tj) to 
just one rid sector: rid — for 4> > 0, or n d — 2 for 
(j) < 0. Then substituting Eq. into the Friedel sum 
rule (fid)o = 25(0)/tt yields 



lim Qi oc ( 

0—0 



;T = 0) 



2sgm? 



aco^TrpoWd). (72) 



The latter relation explains the equality of the exponents 
(3 entering Eqs. ([5^)1 and (|7T|) . It should also be noted 
that Eqs. (Hi]), flBBJl, and (J72) together imply that 



Urn lQ? oc (0;T = O) 

0^0 



lim Txc,im P 



(T). (73) 



At the critical point, the response to a small-but-finitc 
potential 4> obeys another power law, 



Hoc I 



6;A = A c ,T = 0)cx l^l 1 / 5 . 



(74) 



This behavior is exemplified in Fig. [F|] for four different 
values of s. 

Figure [H] shows a logarithmic plot of the static local 
charge susceptibility x Cj i oc (T;c<j = 0) vs temperature T 
for bath exponent s = 0.4 and a number of e-b couplings 
straddling A c . In the quantum-critical regime, the sus- 
ceptibility has the anomalous temperature dependence 



Xc,ioc(T; u = 0) cx T~ x for T, < T < T K , 



(75) 



characterized by a critical exponent x. For T -c T*(A), 
the temperature variation approaches that of one or other 
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12: (Color online) Impurity charge 

<3ioc(A, </>; T — 0) vs e-b coupling A — A c for 
four different values of the bath exponent s. All other 
parameters are as in Fig. [5] As A approaches A c from above, 
lim^Q- Qi oc (A, <f>; T = 0) vanishes (left panel) in a power-law 
fashion (right panel) described by Eq. (|71[) . 
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FIG. 13: (Color online) Impurity charge Qi oc (4>; A = A c , T — 0) 
vs local electric potential |<jf>| for four different values of the 
bath exponent s. All other parameters are as in Fig. [5] The 
dashed lines represent fits to the form of Eq. (|74|l . 



of the stable fixed points. In the Kondo phase, the sus- 
ceptibility is essentially temperature independent, signal- 
ing complete quenching of the impurity, and the zero- 
temperature value diverges on approach to the critical 
coupling as 



Xc,ioc(A; uj = T = 0) oc (A c - for A A" 

In the localized phase, by contrast, 

Ql c (\,4>;T = Q) 



(76) 



Xc,\ oc {T, A;w = 0) = lim 

0— >o 



T 



for A > A c and T«T„ (77) 

indicative of a residual impurity degree of freedom. Pre- 
cisely at the critical e-b coupling, Eq. (jT5[) is obeyed all 
the way down to T = 0. 
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FIG. 14: (Color online) Static local charge susceptibility 
Xc,\oc{T\lo = 0) vs temperature T for s — 0.4, U = —2ed = 
0.1, T = 1.0 (see footnote HI), A = 9, N s = 500, N b = 8, 
and for different values of the e-b coupling A straddling the 
critical value A c — 1.02905. 



TABLE III: Static critical exponents 0, 1/5, x, and 7 defined 
in Eqs. (JTTJ and (ff4)| - (f76) ). respectively, for four different val- 
ues of the bosonic bath exponent s. Parentheses surround the 
estimated nonsystematic error in the last digit. 



s 


P 


1/8 


X 


7 


0.2 


2.0005(3) 


0.6673(1) 


0.1997(2) 


0.997(4) 


0.4 


0.7568(2) 


0.4283(2) 


0.4002(4) 


1.0117(6) 


0.6 


0.3923(1) 


0.2501(7) 


0.600(2) 


1.1805(5) 


0.8 


0.2130(1) 


0.1111(1) 


0.800(2) 


1.703(3) 



Table IIIII lists the numerical values of the critical ex- 
ponents j3, 1/(5, x, and 7, for four different sub-Ohmic 
bath exponents s. For each s, these critical exponents 
are identical within estimated error to those of the spin- 
boson and Ising BFK models. In all cases, we find that 
x = s to within our estimated nonsystematic numerical 
error. We also note that for s < |, the value of 7 lies 
close to its mean- field value of 1. It is conceivable that 
the deviations of 7 from 1 are artifacts of the NRG dis- 
cretization and truncation approximations. 



The exponents in Table IIIII obey the hyperscaling re- 
lations 



S 



1 



1 



2/3 = -x), 7 



(78) 



which are consistent with the ansatz 
' I A — A c | |0| 



F = Tf 



fl/ v ' r p(l+x)/1 



(79) 



for the nonanalytic part of the free energy. Such hy- 
perscaling suggests that the quantum critical point is an 
interacting one^ 
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FIG. 15: (Color online) Imaginary part of the dynamical local 
charge susceptibility Xe'locO^iT — 0) vs frequency to for s = 
0.2, U = -2e d = 0.1, T = 0.5 (see Footnote |56|), A = 9, N a = 
500, Nb = 8, and different e-b couplings A < A c on the Kondo 
side of the critical point, which is located at A c — 0.53008. As 
A — > A^r, x", loc( w i T = 0) follows the quantum critical form 
[Eq. J83j] for T* < < T K , where T K is the Kondo scale of 
the pure-fermionic (A = 0) problem. 
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FIG. 16: (Color online) Critical static and dynamical re- 
sponse: Xc,ioc(T;A = A c ,tj = 0) vs T (circles) and 
Xcioc( a; i A = ^c,T = 0) vs lu (squares) for two representative 
bath exponents s = 0.2 and s = 0.8. All other parameters 
are as in Fig. [5] The equality of the slopes of the static and 
dynamical charge susceptibilities for a given bath exponent 
s indicates that the corresponding critical exponents satisfy 
x = y. 



2. Dynamical local charge susceptibility 
The dynamical local charge susceptibility is 

/>oo 

Xcioc(w,T) = i/ dte- iut (([n d (t)-h M0)-1]». (80) 



Its imaginary part Xcioc can be calculated within the 
NRG as 

Xc,ioc( w ' T ) = Y^r, \( m '\ fl d - l|"i)| 2 



x e 



)S{uj-E m , +E m ). (81) 



Here, |m) is a many-body eigenstate with energy E m , 
and Z(T) = ^ m e~~ Em l T is the partition function. Equa- 
tion (|8"Tj) produces a discrete set of delta- function peaks 
that must be broadened to recover a continuous spec- 
trum. Following standard procedure^ we employ Gaus- 
sian broadening of delta functions on a logarithmic scale: 



S(\o 



\AE\ 



V^b\AE\ 



exp 



(ln|w| -ln|A£|) 2 



b 2 



(82) 

with the choice of the broadening width b = 0.5 In A. 

(a) Zero temperature. Figure [TBI plots XciocC^i^ = 0) 
vs lu for bath exponent s = 0.2 and a series of e-b cou- 
plings A < A c . Whereas X c ',ioc( w ; ^ — — 0) oc to 
for \u>\ -C Tk (the usual Kondo result), we find that 
x"i oc (w; < A < A C) T = 0) oc \uj\ s sgn(uj) as w -> 0, cor- 
responding to a long-time relaxation behavior Xc,ioc(i) oc 
£~(t+ s ). Precisely at the critical e-b coupling, 

Xc,ioc( w ; A = A c , T = 0) ex |wp a sgn(w) for co < T K . 

(83) 



Figure [Tol shows x c 'ioc( CJ :^ = A C ,T = 0) vs w and 
Xc,ioc(r;A — X c ,ui — 0) vs T for representative bosonic 
bath exponents s = 0.2 and s — 0.8. These and all other 
data that we have obtained are consistent with the rela- 
tion 



y 



for < s < 1. 



(84) 



For small deviations from the critical coupling, 
X c ' loc (w; T = 0) exhibits the critical behavior of Eq. (|83|) 
over the range T* <§C |w| <C Tk, where T* is identical (up 
to a constant multiplicative factor) to the crossover scale 
defined in Sec. lIV Cl that vanishes at the quantum critical 
point according to Eq. (|58|) . 

(b) Finite temperatures. Equation (|84| is consistent 
with the presence of lu /T scaling in the dynamical local 
charge susceptibility at the quantum critical point, viz 



X '! loc (u,T;\ = \ c ) =T-** s (u ) /T). 



(85) 



Figure [T71 shows the collapse of data for x" i oc (w, T; A = 
A c ) onto a single function of lo/T within the critical 
regime. The Kondo temperature Tk of the Anderson 
model obtained by setting A = serves as a nonuniversal 
high-frequency cutoff on the critical behavior; the curves 
have a common form for lu/T <C Tk/T. It should be 
noted that the NRG method is unreliable^^ , for \co\<T, 
preventing demonstration of complete u> /T scaling. 

Both the hyperscaling of the static critical exponents 
and what seems to be lo/T scaling of the dynamical 
susceptibility are consistent with the QPT between the 
Kondo and localized phases taking place at an interacting 
critical point below its upper critical dimension. 
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FIG. 17: (Color online) Scaling with uj/T of the imaginary 
part of the dynamical local charge susceptibility ~Xc\ oc (u,T) 
at the critical e-b coupling A c ~ 0.53008 for s = 0.2, U = 
-2e d = 0.1, T = 0.5 (see footnote |g), A = 9, N s = 500, 
N b = 8, and different temperatures T < T K = 0.425. 
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FIG. 18: (Color online) Impurity spectral function A(u;T = 
0) vs frequency u for s = 0.8, U = ~2e d = 0.1, T = 0.01, A = 
3, N s = 1200, Nt — S, and different values of the e-b coupling 
A. For these parameters, C/ e ff defined in Eq. (|33[) changes sign 
at A c o ~ 0.369 and the critical coupling is A c ~ 0.474. 



F. Impurity spectral function 

We now turn to discussion of the impurity spectral 
function j4 ct (o;,T) = — 7r _1 ImGd ff (w, T), where the re- 
tarded impurity Green's function is 

/>oo 

G d<T (u;,T) = -i ( fte <wt «[d ff (t),4(0)] ». (86) 
Jo 

The spectral function can be calculated within the NRG 
using the formulation 



1 



Z(T) 



^2 \(m'\dl\m)\ 



where the notation is the same as in Eq. (|8ip . To recover 
a continuous spectrum, we have again applied Eq. {82]) 
to the delta- function output of Eq. (|57|) . choosing the 
broadening factor b = 0.55 In A that best satisfies the 
Fermi-liquid result A a (u — 0, T — 0) = l/nT for the 
Anderson model. In order to achieve satisfactory results, 
we find it necessary to work with a smaller discretization 
parameter (A = 3 instead of the value A = 9 employed 
for all the quantities reported above) and to retain more 
states (N s = 1200 rather than the 500 that typically 
suffices). Since the spectral functions shown below are 
all spin- independent, we henceforth drop the index a on 
A a . For the particle-hole-symmetric model considered 
in this section, the spectral function is symmetric about 
u = 0. 

Figure [TBI plots A a {u;T = 0) vs oj for s = 0.8 and a 
series of A values. For A = 0, we recover the spectral func- 
tion of the Anderson model, featuring a narrow Kondo 
resonance centered at zero frequency and broad Hubbard 
satellite bands centered around oj = ±if7. Increasing the 
e-b coupling from zero has two initial effects — a displace- 
ment of the Hubbard bands to smaller frequencies, and a 



broadening of the low-energy Kondo resonance — that can 
both be attributed to the boson-induced renormalization 
of the Coulomb interaction described in Eq. ([46]) . 

We expect the Hubbard peak locations to obey ujh — 
±5^cff for < A A c0 . However, the peak loca- 
tions plotted in Fig. [TWa) are better fitted by \ojh\ — 
0AU — \ 2 / (its), which (given the discretization and trun- 
cation effects discussed in Sec. lIII A|) appears to represent 
a stronger bosonic renormalization than that predicted 
by \ujh\ — 7;U c ff. We believe that this discrepancy arises 
primarily from the rapid broadening of the Kondo reso- 
nance with increasing A, which shifts the local maximum 
of the combined spectral function (the sum of the Kondo 
resonance plus Hubbard satellite bands) to a frequency 
smaller in magnitude than the central frequency of the 
Hubbard peak by itself. 

The width 2Fk of the Kondo resonance, plotted in Fig. 
TWb). proves to be equal (up to a multiplicative con- 
stant) to the crossover scale defined in Sec. IIV CI For 
A < A c o, the variation in both scales is well described 
by the replacement of U in the expression 20 for the 
Kondo temperature of the symmetric Anderson model by 
U NRG (U/2) [given by Eq. JUJ], the effective Coulomb 
interaction on entry to the local-moment regime. The 
dashed line in Fig. [T97 b) shows that the resulting for- 
mula, 



r K 



8t/ NRG r 

Ck \ A eX P 



ttA a U 



NRG 



sr 



where Aa is defined in Eq. (|3ip . provides an excellent 
description of Tk over almost the entire range < A < 
A c o — 0.369. This echoes the finding in the Anderson- 
Holstein model that a weak e-b coupling serves primarily 
to reduce the impurity on-site repulsion, leading to an 
increase in the Kondo scaled 

Once the e-b coupling exceeds A c q, further increase in 
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FIG. 19: (Color online) Variation with e-b coupling A of 
two characteristic energy scales extracted from the zero- 
temperature impurity spectral function. All parameters ex- 
cept A are the same as in Fig. 1181 (a) Location ujh of the upper 
Hubbard peak. The dashed line shows u>h = 0.4(7 — A 2 /{ns). 
(b) Kondo resonance width (full width at half height) ZTk- 
The dashed line, representing the prediction of Eq. (188[) 
with C K = 0.82 and with j7 NRG in Eq. ((48} evaluated at 
E — (7/2 — \td.\, fits the data over almost the entire range 
< A < A c0 ~ 0.369. 
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FIG. 20: (Color online) Detail of the impurity spectral func- 
tion A(u>;T = 0) around frequency u = for s = 0.8, 
U = -2e d = 0.1, T = 0.01, A = 3, N 3 = 1600, N b = 8, 
and different e-b couplings A straddling the critical value 
A c ~ 0.47458. For A < A c , A(lu; T = 0) is pinned to the value 
predicted by Fermi-liquid theory. For A > A c , the Kondo reso- 
nance disappears, leaving a pair of low-energy peaks centered 
at \uj\ of order the crossover temperature T* (~1.4 x 10 -8 for 
A = 0.475). 



A leads to suppression of the Hubbard peaks (e.g., see 
the curves for A = 0.4 and A = 0.43 in Fig. [15)1 and to a 
rapid narrowing of the Kondo resonance [see Fig. [TWb)]. 
In the Anderson-Holstein model, the Kondo scale re- 
mains nonzero — although exponentially reduced — for ar- 
bitrarily large e-b couplings. - In the charge-coupled BFA 
model, by contrast, the Kondo peak collapses and Tk ex- 
trapolates to zero as A approaches its critical value A c . 
As shown in Fig. 1201 the central peak remains pinned to 
the Fermi-liquid result A(w = 0, T = 0) = I/71T even as 
the peak width vanishes for A — > \~ . 

In the localized phase (A > A c ), there is no vestige 
of the Kondo resonance, but high-energy Hubbard-like 
peaks reappear; see the curves for A = 0.5 and 0.6 in 
Fig. HU In addition, there is a pair of low-energy peaks 
centered at lo ~ ±T*, as shown in Fig. 1201 



G. Spin-Kondo to charge-Kondo crossover 

Based on the analysis of the zero-hybridization limit 
presented in Sec. IIII A[ one expects spin fluctuations to 
dominate the impurity behavior in the region A -C A c o, 
but charge fluctuations to be dominant for A c o <C A < A c . 
This picture is supported by the behaviors of the ther- 
modynamic susceptibilities discussed in Sec. IIVDI The 
evolution from a spin-Kondo effect to a charge-Kondo 
effect can also be probed by comparing the static local 
charge susceptibility [Eq. ([70]) ] with its spin counterpart 

X s ,ioc(T; uj = 0) = - hm — , (89) 



where h is a local magnetic field that enters an additional 
Hamiltonian term 



■Hsjoc = 2 ("dt -fidi)- (90) 

In particular, characteristic energy scales for the spin and 
charge Kondo effects are expected to be l/x Sl i oc (w = 
0, T = 0) and 4/x c ,ioc(^ = 0, T = 0), respectively [where 
the factor of 4 accounts for the difference in conven- 
tions that tf> couples to fid — 1, whereas h couples to 
(ridf — hdi)/2] ■ Figure I2T1 plots the A dependence of these 
quantities for the parameter set illustrated in Figs. llSl and 
ll~9"l The Kondo resonance width 2Tk crosses over from 
paralleling 1/Xs,loc(0, 0) for small A to loosely tracking 61 
4/Xc,ioc(0, 0) as A approaches A c . In the intermediate 
region near A = A c o, 2Fk is much smaller than either 
inverse static susceptibility, indicating that the Kondo 
effect has mixed spin and charge character. 

Figure [22] presents a A-r phase diagram for s = 0.8 
and fixed U = —2ed, showing data points along the phase 
boundary A = A c (r) and along the crossover boundary 
A = Ax(r), defined as the e-b coupling at which the 
Kondo resonance width 2Fk is maximal for the given 
r. The fact that the latter line rises almost vertically 
from A = A c o at T = provides further confirmation of 
the picture of a crossover from a spin-Kondo effect to 
a charge-Kondo effect resulting from the change in the 
sign of U e g, and establishes the validity of the schematic 
phase diagram (Fig. [3|) presented in the introduction to 
this section. 
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FIG. 21: (Color online) Variation with e-b coupling A < A c of 
the Kondo resonance width 2Fk, the inverse static local spin 
susceptibility l/xs,k>c(w = 0,T = 0), and the inverse static 
local charge susceptibility 4/xc,ioc(w = 0, T = 0). The results 
shown are for s = 0.8, U = -2e d = 0.1, T = 0.01, A = 3, 
N s — 1200, and Nb = 8. For the calculation of the static 
local spin susceptibility via Eq. (|89p . the total spin 5* is not a 
good quantum number, so N s specifies the number of (S Z ,Q) 
states retained after each iteration. 




FIG. 22: Phase boundary A c (r) and crossover boundary 
A x (r) (defined in the text) for s = 0.8, U = -2e d = 0.1, 
A = 3, N s = 1200, and Nt = 8. The data are consistent with 
the schematic phase diagram shown in Fig. [3] 



V. RESULTS: SYMMETRIC MODEL WITH 
OHMIC DISSIPATION 

This section presents results for Hamiltonian ([T]) with 
U = —2£d > and an Ohmic bath (i.e., s = 1). We first 
discuss the behavior of the static local charge suscepti- 
bility. We show that, in contrast with the sub-Ohmic 
case < s < 1, the crossover scale vanishes in exponen- 
tial (rather than power-law) fashion as the e-b coupling 
approaches its critical value from below, and there is no 
small energy scale observed on the localized side of the 
transition. Therefore, the QPT for the Ohmic case is 
of Kosterlitz-Thouless type. At the end of the section, 
we study the effects of the e-b coupling on the impurity 




FIG. 23: Schematic renormalization-group flows on the A-A 
plane for the symmetric model with bath exponent s = 1. 
Trajectories represent the flow of the couplings A entering 
Eq. (|15p and A defined in Eq. (|50[) under decrease in the 
high-energy cutoffs on the conduction band and the bosonic 
bath. A separatrix (dashed line) forms the boundary between 
the basins of attraction of the Kondo fixed point (K) and a 
line of localized fixed points (L). Flow along the separatrix is 
toward the free-orbital fixed point (FO) located at A = A c o. 
For A = only, there is flow away from FO toward the local- 
moment fixed point (LM) at A = 0. 



spectral function. 



A. Fixed points and thermodynamic 
susceptibilities 

Figure [23] plots the schematic renormalization-group 
flows for a symmetric impurity coupled to an Ohmic 
bath. The flows within the Kondo basin of attraction 
are qualitatively very similar to those for the sub-Ohmic 
case depicted in Fig.UJ In the localized regime, however, 
the e-b coupling flows not to A = oo, but rather to a fi- 
nite limiting value that varies continuously with the bare 
values of A and T. What is shown as a line of fixed points 
in Fig. [23] is really a plane of fixed points described by 
Hlc [Eq. ([49]) ] with effective couplings A > A c0 , W p = 0, 
and < Wd < oo. Another important departure from 
the sub-Ohmic case is that for s = 1 there is no longer 
a distinct critical point reached by flow along the sepa- 
ratrix from the free-orbital fixed point; rather these two 
fixed points merge as s — > 1 — , leaving a critical end- 
point at A = A c o, A = 0. Strictly, this is a line of critical 
endpoints described by H^c [Eq. P5]) ] with effective cou- 
plings A = A c o, W p = 0, and < Wd < oo. For a fixed 
bare value of T, the endpoint value of Wd is just the 
limit of the localized fixed-point value of Wd as the bare 
coupling A approaches the phase boundary A c (r). 

The behaviors of the static impurity spin and charge 
susceptibilities are qualitatively very similar to those for 
a sub-Ohmic bath, as discussed in Sec. IIVDI The only 
significant difference is that for s = 1, limj^o Tx c ,imp(T) 
undergoes a discontinuous jump from its value of for 
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FIG. 24: (Color online) Static local charge susceptibility 
Xc,1oc{T;lu = 0) vs temperature T for s = 1, U = —2e d — 0.1, 
r = 0.01, A = 9, N s = 800, N b = 12, and different e-6 cou- 
plings A. On the Kondo side of the QPT (A < A c ~ 0.726), 
there is a clear crossover from quantum-critical to screened 
behavior around the renormalized Kondo temperature T* = 
4/xc,ioc(^ = T — 0). No such crossover is evident on the 
localized side (A > A c ). 

A < A c to a nonzero value for A = A+. This jump can be 
understood through Eqs. ([tPo]) and (prj) as a consequence 
of the fact that Wd does not diverge on approach to the 
critical coupling. 

B. Static local charge susceptibility and crossover 
scale 
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FIG. 25: Variation with e-b coupling A of the local charge 
susceptibility Xc,loc(w = T = 0) in the Kondo phase A < A c ~ 
0.726 and of the order parameter lim^^ - Qioc(<^;T = 0) in 
the localized phase A > A c , for s = 1, U = —2ed — 0.1, 
r = 0.01, A = 9, N s = 800, and N b = 12. The dotted line 
shows a fit of the susceptibility data using Eqs. ([91) and (|92j. 



magnitude of this jump is nonuniversal, being related via 
Eq. ([72]) to the value of Wd at the critical endpoint. 

The properties described above are analogous to those 
of the Kondo model [Eq. (J2U)] at the transition between 
the Kondo-screened phase (reached for J± ^ and J z > 

— \J±\) and the local- moment phase (reached for J z < 

— | Jj_|). Such behaviors are characteristic of a Kosterlitz- 
Thouless type of QPT. 



Figure[24lis a logarithmic plot of the static local charge 
susceptibility Xc,\oc(T;lo = 0) vs temperature T for dif- 
ferent e-b couplings A. On the Kondo side of the phase 
boundary, Xc}oc(T;lu — 0) is proportional to \/T at high 
temperatures, but levels off for T < T*. We find it con- 
venient to define 



T* = 4/xc,k>c(u> = T = 0) for A -> X~ 



(91) 



thereby removing the ambiguity in the definition of the 
crossover iteration (see Sec. IIV C| on the Kondo side 
of the s = 1 quantum phase transition. 

For A — » A^ , the crossover scale vanishes according to 
(see Fig. [25) 



T* oc exp 



V 7 ! - (VAc) 2 



(92) 



In the localized phase, Xc,\oc{T\lo = 0) satisfies Eq. 
(|77jl over the entire temperature range T <C U . Since the 
critical and localized fixed points share the same temper- 
ature variation, no crossover scale can be identified on 
the localized side of the phase boundary. Moreover, the 
order parameter lim0_» o Qioc(</ ) ; T = 0) does not vanish 
continuously as A — ► A+, but rather undergoes a discon- 
tinuous jump at the transition, as shown in Fig. 1251 The 



C. Impurity spectral function 

Figure [21] shows the impurity spectral function 
A(ui;T = 0) for an Ohmic bath. The behavior in the 
Kondo phase is similar to that in the sub-Ohmic case 
discussed in Sec. IIV Fl As the e-b coupling A increases 
from zero, the Hubbard satellite bands are initially dis- 
placed to smaller frequencies according to u>h — ±\U c s 
[Fig. [2"7f a)]. while the width 2Fk of the Kondo resonance 
[Fig. l2"7T b)] first rises before falling sharply on approach 
to A = A c . Just as for < s < 1, the variation in Tk 
for A < A c0 is well described by Eq. ([55]) with L> NRG 
[Eq. J48J] evaluated at E = U/2. Throughout the Kondo 
phase, A{oj = T = 0) remains pinned at its Fermi-liquid 
value I/'kY. 

For A > A c , however, the behavior of the spectral func- 
tion is quite different for s = 1 than for < s < 1. In 
the sub-Ohmic case, the Kondo-phase pinning extends 
to the quantum critical point, i.e., t:YA(llj = T = 0, A = 
A c ) = 1, while in the localized phase peaks appear at 
uj ~ ±T». Figure [28l shows that the Ohmic spectral func- 
tion instead satisfies TrYA{u> = T = 0, A = A c ) < 1, and 
exhibits no feature in the localized phase at energy scales 
much smaller than ^|[/ c ff|- 
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FIG. 26: (Color online) Impurity spectral function A(ui;T — 
0) vs (j for s = 1, U = -2e d = 0.1, T = 0.01, A = 3, N s = 
1200, Nb = 12, and different values of the e-b coupling A. For 
these parameters, U e s [Eq. (|33p] changes sign at A c o — 0.413 
and the critical coupling is A c ~ 0.669. 



0.016 



0.040 o-e 




0.004 



0.030 



0.2 0.4 0.6 
X 



FIG. 27: (Color online) Variation with e-b coupling A of 
two characteristic energy scales extracted from the zero- 
temperature impurity spectral function. All parameters ex- 
cept A are the same as in Fig.(5S] (a) Location ljh of the upper 
Hubbard peak. The dashed line shows aJfr(A) = 0AU — A 2 /it. 
(b) Kondo resonance width (full width at half height) 2Fk. 
The dashed line, representing the prediction of Eq. 1)880 
with C K = 0.82 and with [7 NRG in Eq. (gSJ evaluated at 
E — U/2 — |e<j|, fits the data over almost the entire range 
< A < Aco ^ 0.413. 



VI. RESULTS: ASYMMETRIC MODEL 

Sections IIVI and [V] focused exclusively on results for 
a symmetric impurity satisfying td = —U/2 in Eq. @ 
or, equivalently, 5d — in Eq. ©. We now turn to 
the general situation of an asymmetric impurity, starting 
with the sub-Ohmic case < s < 1. 

For 5d 7^ and small, nonzero values of A, one expects 
the fermionic sector of the charge-coupled BFA model to 
behave in essentially the same manner as in the asym- 
metric Anderson model (reviewed in Sec. lIIIBj) . with the 
exception that the effective value of the Coulomb inter- 
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FIG. 28: (Color online) Impurity spectral function A(u;T = 
0) vs frequency u) on a logarithmic scale for s = 1, U = 
-2e d = 0.1, T = 0.01, A = 3, N s = 1200, N b = 12, and 
different e-b couplings A. For A < A c ~ 0.669, the behavior is 
similar to that found for < s < 1. However, for A > A c , the 
spectral function is essentially featureless below the energy 
scale ^|[/ e ff| of the Hubbard peaks. 



action U will be reduced by the coupling to the bosonic 
bath. At temperatures well below 7k, there will be no 
further renormalization of the electronic degrees of free- 
dom, the system will exhibit quasiparticle excitations de- 
scribed by Hg£ G in Eq. and the low-energy many- 
body states will share a nonvanishing expectation value 
(fid — 1) [= Qioc{T = 0)]. The bosons will couple to 
this impurity charge, yielding low-energy states described 
most naturally in terms of displaced-oscillator states [cf. 
Eq. jZrj] annihilated by operators 




h d - 1) 



(93) 



For s < 1, the e-b coupling is relevant so A will scale to 
strong coupling below a crossover temperature Tl <Si Tk- 

For 5d 7^ and very large values of A, one instead ex- 
pects the bosons to localize the impurity at a high tem- 
perature scale Tl into a state with (rid) ~ (for Sd > 0) 
or (rid) ~ 2 (for 5 d < 0). For T < T L , the impurity de- 
grees of freedom will be frozen, the bosonic spectrum will 
rapidly approach strong coupling, and the conduction 
electrons will have an excitation spectrum correspond- 
ing to H^ G in Eq. ([35]) with a small value of \V \. 

Given the equivalence of and i?pj RG , it seems 

likely that the low-energy behavior of the asymmetric 
model will be the same in the small-A and large-A limits. 
This suggests that the many-body eigenstates evolve adi- 
abatically as the e-b coupling is increased from A = + to 
A — > oo, without the occurrence of an intervening QPT. 

For s = 1, the e-b coupling is marginal, rather than rel- 
evant. One again expects a continuous evolution of the 
low-energy NRG spectrum with the bare value of A. How- 
ever, in this Ohmic case, the bosonic excitations should 
correspond to nonintcracting displaced oscillators rather 
than the (truncated) strong-coupling spectrum found for 
< s < 1. 
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FIG. 29: (Color online) Variation in the magnitude (1 — hd)o 
of the ground-state impurity charge with e-b coupling A for 
s = 0.4, U = 0.1, T = 0.01, A = 9, N s = 500, and N b = 8. 
Symbols represent results for five values of the impurity asym- 
metry Sd — £d + U/2. The solid lines corresponding to each 
case Sd ^ represent the impurity charge calculated by solv- 
ing the Anderson model [Eq. (II ip ] for the same Sd value but 
using an effective Coulomb interaction [7 NRG (0.3t/) [Eq. (gHJ)]- 
The Sd = symbols show values of lining- <3ioc(A, (j>;T = 0), 
connected by an interpolating line. 

The preceding arguments are supported by our NRG 
results. Here, we illustrate just the case s — 0.4. Fig- 
ure [29] shows the variation with A of the ground-state 
expectation value (1 — fid)o for several values of 8 d . In 
the symmetric case (the Sd = curve in Fig. |2"9"]) . the 
impurity charge vanishes throughout the Kondo phase, 
and grows in power-law fashion on entry to the localized 
phase. Away from particle-hole symmetry, by contrast, 
(1 — hd)o increases smoothly from its Anderson-model 
value at A = to approach 1 as A — * oo. 

For all nonzero values of Sd, T and A, the low-energy 
spectrum can be decomposed into the direct product of 
the fermionic spectrum corresponding to HgQ G (Vi) [or 
Hpi (Vb)] and the same localized-phase bosonic spec- 
trum as found for the symmetric model. The potential 
scattering V\ (or Vq) is tied to {fid — l)o by Eq. ([3"T]) , just 
as in the Anderson model. 

For small A, the value of (fid — l)o can be related to 
the corresponding quantity in the Anderson model by 
making use of the effective Coulomb interaction intro- 
duced in Sec. IIII Al In the asymmetric Anderson model, 
the ground-state charge becomes frozen once the sys- 
tem passes out of its mixed- valence regime, i.e., some- 
what below a characteristic temperature Tf defined^ for 
r <C — £rf <C U as the solution of 

T f = \e d \- -In ^. (94) 

7T If 

In the charge-coupled BFA model, U and in Eq. 
(IM|) should presumably be replaced by U(Tf) and Sd — 
hU(Tf), respectively. However, it suffices for our pur- 
poses to note that Tf can be expected to be of the same 
order as, but somewhat smaller than, |e<j|. It is then rea- 
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FIG. 30: (Color online) Variation in the bosonic localiza- 
tion temperature Tl with coupling A for s = 0.4, U = 0.1, 
T = 0.01, A = 9, N s = 500, Nb = 8, and various impurity 
asymmetries Sd — £d + U /2. The solid lines were obtained by 
evaluating Eq. fjQ5p with the (1 — hd)o values shown in Fig. 
[29] and with C L = 3. 



sonable to hypothesize that (rid ~ l)o i n the asymmetric 
charge-coupled BFA model should be close to the ground- 
state impurity charge of the Anderson model with the 
same T and S d , but with U replaced by U(E) [Eq. ([47]) ] 
evaluated at E ~ Tf. Our numerical results support 
this conjecture. For example, Fig. [29] shows that close 
to particle-hole symmetry (e d = ~U/2), the Anderson- 
model charge calculated for U NRG (E) [Eq. ggj] with 
E = 0.3C/ (solid lines) reproduces quite well the value 
of (fid — l)o (symbols) over quite a broad range of e-b 
couplings < A < |A C , where A c ~ 0.29835 is the criti- 
cal coupling of the symmetric problem. 

In the small- A limit, one can also estimate the boson- 
localization temperature Tl by considering the evolution 
with decreasing T of the effective value of X(hd — l)o- The 
impurity charge does not renormalize, while to lowest or- 
der the effective e-b coupling obeys^ Eq. ([28]). Defining 
Tl by the condition A(T£,)|(n<j — 1) o | = Cx, we find 

TL^\CZ l \{h d -l) \ V{1 ~ s) . (95) 

In Fig. 1301 symbols represent Tl values extracted from 
the crossover of bosonic excitations in the NRG spec- 
trum, while solid lines show the results of evaluating Eq. 
(|95[) using Cl = 3 and the (fid — l)o values shown in 
Fig. [211 The algebraic relation between the numerical 
values of Tl and (1 — n d )o is well obeyed over a range of 
e-b couplings that extends beyond A c of the symmetric 
problem. 

Figure [31] plots the static local charge susceptibility 
calculated for s = 0.4 at the critical e-b coupling of the 
symmetric model. For Sd ^ 0, Xc.ioc follows the quan- 
tum critical behavior Xc,\oc(T] u) — 0) cx T~ x from a 
high-temperature cutoff of order Tk down to a crossover 
temperature T*, below which the susceptibility saturates. 
Based on Eq. ([79"]) with the identification cj>= Sd, one ex- 
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FIG. 31: (Color online) Static local charge susceptibility 
Xc,ioc(-T;<^ = 0) vs temperature T for s = 0.4, U = 0.1, 
r = 0.01, A ~ 0.29835, A = 9, N s = 500, N b = 8, and 
various impurity asymmetries Sd = ed + U/2. The e-b cou- 
pling equals the critical coupling A c of the symmetric case 
Sd = 0. Inset: zero-temperature static local charge suscepti- 
bility Xc,ioc(w = T = 0) vs 6d- 



pects T* oc |^d| 2 ^ 1+:r ' and, hence, 

XciocW>; A = A c , w = T = 0) oc \5 d \- 2x /( 1+x l (96) 

The log-log plot in the inset of Fig. [31] has a slope 0.57 
that is fully consistent with Eq. (|96|) . 

The results of this work show that gaining direct ac- 
cess to the quantum critical point of the charge-coupled 
BFA model requires simultaneous fine tuning of two pa- 
rameters: the e-b coupling A as a function of the hy- 
bridization r and the on-site Coulomb repulsion U\ and 
the particle-hole asymmetry (determined in our calcula- 
tions solely by 8 d — e d + U/2, but in general also affected 
by the shape of the conduction-band density of states). 
While it may prove very challenging, or even impossible, 
to achieve this feat in any experimental realization of the 
model, it should be a more feasible task to carry out a 
rough tuning of parameters that places the system in the 
quantum critical regime over some window of elevated 
temperatures and/or frequencies. 



duction band and is coupled, via its charge, to a dissipa- 
tive environment represented by a bosonic bath having 
a spectral function that vanishes as ui s for vanishing fre- 
quencies uj — > 0. With increasing coupling between the 
impurity and the bath, we find a crossover from a conven- 
tional Kondo effect — involving conduction-band screen- 
ing of the impurity spin degree of freedom — to a charge- 
Kondo regime in which the delocalized electrons quench 
impurity charge fluctuations. 

Under conditions of strict particle-hole symmetry, fur- 
ther increase in the impurity-bath coupling gives rise for 
< s < 1 to a quantum phase transition between the 
Kondo phase, in which the static charge and spin sus- 
ceptibilities approach constant values at low tempera- 
tures, and a localized phase in which the static charge 
susceptibility exhibits a Curie- Weiss behavior indicative 
of an unquenched local charge degree of freedom. For 
sub-Ohmic bosonic bath spectra (described by an expo- 
nent s satisfying < s < 1), the continuous quantum 
phase transition is governed by an interacting critical 
point characterized by hyperscaling relations of critical 
exponents and uj/T scaling in the dynamical local charge 
susceptibility. Moreover, the continuous phase transition 
of the present model belongs to the same universality 
class as the transitions of the spin-boson and the Ising- 
anisotropic Bose-Fermi Kondo models. For an Ohmic 
(,s = 1) bosonic bath spectrum, the quantum phase tran- 
sition is of Kosterlitz-Thouless type. 

In the presence of particle-hole asymmetry, the quan- 
tum phase transition described in the previous para- 
graph is replaced by a smooth crossover, but for small- 
to-moderate asymmetries, signatures of the symmetric 
quantum critical point remain in the physical proper- 
ties at elevated temperatures and/or frequencies. Inves- 
tigation of the regime of strong particle-hole asymme- 
try, and of self-consistent versions of the charge-coupled 
Bose-Fermi Anderson model that arise with the extended 
dynamical mean-field theory of lattice fermions, will be 
pursued in future work. 
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